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INTRODUCTION 

The  problem  of  providing  safe  enclosures  for  the  public  from  fallout 
radiation  resulting  from  a  nuclear  explosion  involves  the  analysis  of  a 
variety  of  structures  to  determine  what  protection  each  affords  from  the 
radiation.   This  analysis  may  be  accomplished  either  theoretically  or  ex- 
perimentally.  A  theoretical  method,  the  "Engineering  Manual"  method,  has 
been  tested  extensively  and  has  been  shown  to  give  good  results  in  most  cases. 
The  experimental  methods  which  simulate  a  fallout  field  around  a  full- sized 
structure  yield  better  results,  but  are  cumbersome,  expensive,  and  time- 
consuming.   Because  of  the  need  to  more  fully  check  the  theoretical  tech- 
niques, and  because  of  the  disadvantages  of  using  the  experimental  method 
on  the  full-scale  buildings,  scale  models  of  the  structures  of  interest  are 
being  used. 

Modeling  techniques  have  been  developed  to  a  large  extent  and  have 
proven  to  be  useful.   However,  much  work  remains  to  be  done  before  model 
results  can  be  accepted  as  valid  in  all  cases  on  their  own  merits.   In  this 
work  new  and  original  model  data  were  collected  in  an  effort  to  help  ac- 
complish this  goal.   A  comparison  has  been  made  of  protection  in  1:12  and 
l:k   scale  models  of  a  concrete  block  house.   The  results  have  been  compared 
with  theoretical  predictions  and  with  previous  studies  made  in  the  full-rcale 
structure. 


THEORY  OF  MODELING 

The  radiation  dose  distribution  inside  a  structure  due  to  a  radioactive 
field  outside  in  theory  should  be   exactly  reproduced  in  a  scale  model  of  the 
building  provided  the  following  two  requirements  are  met: 

1.  The  density  of  all  materials  affecting  the  distribution  of 
radiation  must  be  increased  by  the  same  factor  by  which  all 
dimensions  are  reduced. 

2.  The  gamma-ray  scattering  and  absorption  properties  per  unit 
mass  of  the  materials  used  in  the  scaled  experiment  must  be 
the  same  as  those  in  the  full-scale  experiment. 

These  conditions  require  that  the  densities  of  the  air  and  ground  as  well  as  the 
building  materials  be  increased  by  the  same  factor  that  linear  dimensions  are 
reduced.   In  practice  the  density  of  these  materials  can  not  be  increased  by  a 
factor  of  10,  which  is  necessary  if  the  scale  model  is  to  provide  much  advantage 
over  the  full-size  structure. 

Since  it  is  impractical  to  satisfy  the  above  requirements,  approximations 
must  be  made.   The  first  is  to  replace  the  building  material  with  one  more 
dense.   For  this  work  steel  was  used  to  replace  concrete  since  it  is  high 
density,  inexpensive,  and  not  too  different  in  nuclear  properties.   This  per- 
mitted an  increase  in  density  about  a  factor  of  four.   Prior  experiments  have 
shown  that  scaling  is  still  realistic  if  wall  thicknesses  are  not  greater  than 
10  percent  of  the  average  dimensions  of  a  given  room.   This  allowed  the  total 
mass  of  the  walls  to  be  increased  by  another  factor  of  three.   Consequently,  a 
scaling  factor  of  12  was  achieved. 

There  are  three  possible  criteria  for  selection  of  model  wall  thicknesses: 


1.  The  mass  thickness  of  the  walls  of  the  full-size  structure 
may  be  duplicated. 

2.  The  electron  density  may  be  maintained. 

3.  The  broad-beam  absorption  data  for  slabs  may  be  used. 

To  illustrate  the  difference,  calculations  have  shown  that  1.69"  of  iron  is 
necessary  to  maintain  the  same  mass  thickness  as  an  8"  thick,  light-weight 
concrete  wall;  1.80"  is  required  if  the  electron  densities  are  matched;  and 
1.6h"  is  needed  if  the  broad-beam  data  are  used. 

Since  the  use  of  the  broad-beam  data  places  a  heavy  dependence  on  the 
linearity  of  the  detectors  with  gamma-ray  energy,  criterion  3  is  seldom  chosen 
for  determination  of  barrier  thicknesses. 

The  average  cobalt-60  gamma-ray  energy,  1.25  Mev. ,  is  representative  of  the 
fallout  radiation  energy  spectrum  at  1.12  hours  after  a  fission  event  and  is 

conservative,  insofar  as  shielding  calculations  are  concerned,  for  fallout  energy 

2  3 
spectra  at  other  times.  '   At  this  energy  Compton  scattering  is  the  most  prob- 
able method  of  interaction  of  gamma  radiation  with  a  shield.   This  effect  de- 
pends upon  the  electron  density  of  the  material.   If  the  barriers  of  the 
structure  under  study  were  thin,  then  the  model  should  be  constructed  so  as  to 
match  the  electron  density  of  the  walls  of  the  full-size  building.   However,  the 
walls  of  the  concrete  block  house  for  this  work  were  69  psf,  which  is  approxi- 
mately two  mean  free  paths  of  cobalt-60  gamma  rays.   Thus,  the  greater  portion  of 
the  radiation  reaching  the  interior  of  the  block  house  has  been  scattered  in  the 
barriers.   The  attenuation  of  the  scattered  radiation  is  more  dependent  on  the 
effective  atomic  number  of  the  shields.  Thus,  the  model  barrier  thicknesses  for 

this  work  were  calculated  on  the  basis  of  maintaj     the  same  mass  thickness 

k 
present  in  the  full-scale  structure.   Experimental  data  suppor' 


While  iron  has  a  smaller  electron  density  than  concrete  and  will  transmit 
more  unscattered  radiation  than  the  same  mass  thickness  of  concrete,  this 
increase  is  offset  by  the  greater  probability  of  absorption  of  scattered 
radiation  in  iron  compared  to  that  of  concrete. 


EXPERIMENTAL  MEASUREMENTS 

Description  of  Models 

One-Fourth  Scale  Model 

The  steel  l:k   scale  model  was  a  replica  of  a  single  story  concrete  block 
house  with  a  basement  (Fig.  11)  located  at  the  Kansas  State  University  Nuclear 
Engineering  Shielding  Facility  (KSUNESF).   The  external  dimensions  of  the 
model  were  kept  very  close  to  l/k   of  the  full-scale  values,  while  wall,  roof, 
and  floor  thicknesses  were  calculated  for  the  model  to  give  the  same  mass 
thickness  as  those  of  the  concrete  block  house.  These  values  were  I.69"  of 
steel  for  the  walls  and  1.35"  for  the  roof  and  floor.   Since  slabs  of  steel 
with  these  thicknesses  were  unavailable  commercially,  values  of  1.75"  and  1.375 
were  used.   These  compare  to  2.0"  of  iron  for  the  walls  and  1.5"  for  the  roof 
and  floor,  had  the  dimensions  of  the  block  house  simply  been  scaled  down  by  a 
factor  of  four. 

Each  wall  consisted  of  four  slabs  of  steel  (Fng.  2).   The  floor  was  three 
layers  thick  with  each  layer  made  up  of  four  slabs  arranged  so  as  to  leave  an 
opening  into  the  basement.   A  trap  door  was  constructed  to  fill  the  entrance. 
Three  layers  of  steel  bolted  together  made  up  the  roof.   Each  layer  consisted  of 
two  slabs  laid  side  by  side  such  that  cracks  in  alternate  layers  were  perpendi- 
cular.  The  three  layers  were  bolted  together  with  four  l/2"  eye  bolts,  and  two 
cables  were  attached  through  the  eyelets.  This  facilitated  the  removal  of  I 
roof  with  a  fork  lift  for  positioning  or  retrieving  dosimeters.   The 
up  the  mode]     !  ,  roof,  and  floor  were  designed  and  assembled  so  th 
could  be  no  streaming  of  gamma-rays  through  cracks  or  o'  ;r  openings,   'i 
was  held  to      by  two  3/^"  steel  bands  strapped  around  i  . 
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Figure  1.   Drawing  of  the  l:k   model  of  the  KSUNESF  block  house. 


Plan  View 


Figure  2.   Detail  of  one  corner  of 
the  l:k   model. 
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A  unique  feature  of  this  model  was  that  it  could  be  adapted  to  simulate  a 
structure  having  a  portion  of  its  "basement  walls  exposed.   This  was  accomplished 
by  jacking-up  the  entire  structure,  adding  the  exposed  "basement  walls,"  and 
then  lowering  the  model  down  onto  the  walls.   Two  steel  rods  held  the  four  base- 
ment walls  rigid  (Fig.  3). 

One-Twelfth  Scale  Model 

The  1:12  scale  model  was  also  constructed  of  steel  and  was  a  replica  of  the 
KSUNESF  block  house.   The  model  dimensions  were  l/l2  the  corresponding  dimensions 
of  the  full-size  structure,  except  for  the  thicknesses  of  the  walls,  roof,  and 
floor  which  were  calculated  for  the  model  to  give  the  same  mass  thicknesses  as 
those  in  the  block  house.   It  was  necessary  that  l/8"  thick  sheets  of  steel  be 
added  to  the  roof,  walls,  and  floor  of  the  original  model  in  order  that  the  mass 
thicknesses  be  closely  matched.   This  made  the  walls  1.75"  and  the  floor  and  roof 
1.375"  thick.   See  Fig.  k   .   As  illustrated  in  Fig.  h, ,   this  model  could  also  be 
adapted  to  simulate  a  building  with  a  portion  of  its  basement  walls  exposed. 

Experimental  Procedure 

The  experimental  procedure,  much  the  same  for  both  models,  was  as  follows: 

1.  Polyethylene  tubing  was  laid  around  the  model  in  a  prescribed  fashion 
to  simulate  a  limited  plane  of  fallout  radiation.   Figs.  5  and  8 
show  the  tubing  layouts  for  the  1:U  and  1:12  models,  respectively. 

2.  The  source  container  was  positioned  behind  the  concrete  block  house 
about  50'  from  the  1:12  model  and  80'  from  the  larger  model. 

3.  The  pumping  system  was  made  ready  and  one  or  more  dummy  source  runs 
were  made. 

h.        Charged  dosimeters  were  enclosed  in  plastic  bags  and  placed  at 

planned  locations  in  the  first  floor  and  basement  of  the  model  being 
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used.   This  was  done  during  the  assembly  of  the  1:12  model  and 
involved  the  removal  of  the  roof  and  trap  door  of  the  1:U  model. 

5.  A  source  with  an  approximate  strength  of  76  curies  was  pumped 
through  a  section  of  tubing  at  a  constant  speed,  so  that,  since 
the  tubing  density  was  constant  in  that  section,  a  uniformly 
distributed  fallout  field  was  simulated.   Details  of  the  operation 
of  the  pumped  source  system  including  safety  measures  are  given  in 
Reference  5. 

6.  The  dosimeters  were  retrieved,  unwrapped,  read,  charged,  rewrapped 
in  plastic  bags  and  returned  to  a  position  in  the  model,  if  another 
run  were  to  be  made. 

7.  Along  with  the  dosimeter  readings,  temperature,  atmospheric  pressure, 
and  the  exposure  time  were  recorded. 

Measurements  and  Results 

One-Fourth  Scale  Model 

The  first  series  of  experimental  measurements  on  the  1:U  model  were  made 
with  the  floor  of  the  model  flush  with  the  ground  level  and  with  the  inner 
section  of  tubing  (Area  II)  laid  out  as  shown  in  Fig.  5.   Dosimeters  were 
placed  in  each  corner,  on  the  sides,  and  in  the  center  of  both  floors  of  the 
model  and  were  positioned  nine  inches  above  the  floor  level  -  the  equivalent 
of  three  feet  in  the  full  scale  structure.   See  Fig.  6.   The  2r  dosimeters  were 
used  in  the  first  floor  and  10  mr  dosimeters  in  the  basement.   After  several 
preliminary  runs  had  been  made,  it  was  found  necessary  to  bake  out  the  10  mr 
chambers  in  an  effort  to  obtain  less  erratic  readings.   The  experiment  was 
then  redone  with  the  baked  out  chambers  in  the  basement.   This  time  more  con- 
sistent results  were  observed  as  shown  by  the  raw  data  in  Table  F-5. 
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Area  IV 

31^2  Sq.  Ft. 

Tubing  interval  =  2 ' 


Area  III 

1368  Sq.  Ft. 

Tubing  Interval  =  1' 


Area  II 

1105  Sq.  Ft. 

Tubing  Interval  =  1' 


Area  I   (See  Fig.  7 )   ~ 


1:U  Model 
5 '  Square 


E  i  fure  5.  TubJ       :  for  l:h   model. 
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Figure  6.   Drawing  showing  the  detector  locations  with  the  basement  submerged 
(A)  and  with  the  basement  partially  exposed  (B).- 
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The  first  run  with  tubing  Area  III  (27* -bQ')    and  the  model  floor  flush 
with  the  ground  level  was  done  with  2r  dosimeters  in  the  first  floor  and  10  mr 
dosimeters  in  the  basement.   Such  low  readings  were  obtained  with  the  2r  dosi- 
meters that  they  were  replaced  by  200  mr  chambers  in  subsequent  runs.   These 
data  are  reproduced  in  Table  F-6. 

For  the  last  tubing  area  and  the  same  model  configuration,  data  could  only 
be  taken  on  the  first  floor  (Tab]  e  F-8)  since  it  would  have  required  an  excessive 
exposure  time  to  get  satisfactory  results  in  the  basement. 

The  contamination  in  Area  I  was  simulated  by  a  series  of  point  sources.   As 
shown  in  Fig.  7,  a  U5   sector  out  to  a  radius  of  five  feet  was  divided  into 
seven  smaller  areas.   A  point  source  was  placed  on  the  centroid  of  each  area  in 
turn  until  all  seven  areas  had  been  covered.   The  exposure  time  for  each  area 
was  calculated  so  as  to  maintain  a  uniform  average  source  density. 

As  might  be  expected,  a  considerable  amount  of  radiation  reached  the 
dosimeters  in  the  basement  directly  by  penetrating  the  concrete  and  earth  adja- 
cent to  the  model.   Since  no  provisions  were  made  for  this  ground  penetration  in 
the  theoretical  calculations  ,  and  since  the  contribution  was  comparatively  small 
in  the  full-scale  block  house,  it  was  necessary  to  eliminate  this  contribution. 
This  was  done  by  performing  the  experiment  with  the  point  source  just  described 
except  that  now  the  point  source  was  covered  with  a  minimum  of  four  inches  of 
lead.  This  eliminated  practically  all  radiation  in  the  basement  except  that 
which  penetrated  the  concrete  and  ground.   The  ground  contribution  from  the  field 
farther  out  than  five  feet  from  the  center  of  the  model  was  negligible. 

With  the  conclusion  of  the  first  series  of  experiments,  the  model  was 
elevated  six  inches  so  as  to  simulate  a  block  house  with  two  feet  of  its  ba 
walls  exposed,  and  the  above  experiments  wer        i.   Even  though  the  first 
floor  of  the  model  was  elevated  six  inches  and  nothi  led  in  the  floor 
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Figure  7.   Area  I  of  the  l:k   Model  Field 
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of  the  basement  to  take  up  the  additional  six  inches  created,  the  basement  was 
considered  to  be  the  same  size.   Thus,  the  dosimeters  in  the  basement  were 
raised  an  additional  six  inches  to  maintain  the  three-foot  level  in  the  full- 
scale  building. 

As  was  the  case  with  the  first  series  of  experiments,  the  source  exposure 
times  were  adjusted  insofar  as  possible  to  maintain  dosimeter  readings  between 
20  and  80  percent  of  their  full-scale  values.   The  raw  data  collected  from 
experiments  on  the  l:k   model  begin  on  Page  59  of  Appendix  F.   The  reduced  data 
and  protection  factors  calculated  therefrom  appear  in  Tables  1  and  3  > 
respectively. 

One-Twelfth  Scale  Model 

The  tubing  for  the  1:12  model  was  laid  completely  around  the  model  and  in 
two  sections,  each  representing  an  annular  strip  of  radioactive  fallout.   The 
field  closest  to  the  model,  Area  I,  was  represented  by  a  series  of  point  sources 
over  a  sector  including  l/8  of  the  perimeter  of  the  model.   This  area  was  sub- 
divided into  seven  smaller  areas  as  shown  in  Fig.  9  ,   ^nd   the  experiment  was 
conducted  as  outlined  for  the  l:k   model.   Since  the  1:12  model  was  easily 
assembled  and  disassembled,  each  series  of  experiments  were  performed  first  with 
the  model  floor  flush  with  the  ground  level  and  then  with  the  floor  elevated. 

Only  the  200  mr  and  2r  dosimeters  were  used  with  the  1:12  model  since  they 
were  the  smallest  physical  size  available.   The  200  mr  dosimeters  were  generally 
used  in  the  basement  while  the  2r  dosimeters  were  required  for  the  first  floor 
measurements.   Both  types  were  positioned  in  the  same  relative  locations  as  in 
the  1:U  model;  that  is,  three  inches  from  the  floors  in  the  corners,  sides,  and 
center  of  the  model. 


IT 


Figure  8 .   Tubing  layout  for  the  1:  12  mod 
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Some  difficulty  was  experienced  collecting  data  from  Area  III  with  the 
basement  of  the  model  completely  underground.   As  can  be  seen  from  Table  F-26, 
lower  than  desirable  doses  were  accumulated  in  the  basement  in  spite  of  the 
lengthy  exposure  times  which  averaged  about  three  hours  for  each  run. 

Since  these  and  several  other  runs  were  rather  long,  all  three  types  of 
dosimeters  were  checked  to  determine  the  extent  of  leakage.   It  was  found  that 
the  200  mr  and  2r  dosimeters  did  not  drift  significantly  in  2k   hours.   The  10  mr 
chambers,  although  not  as  stable,  did  not  drift  enough  to  require  corrections  to 
the  experimental  data. 

As  with  the  l:k   model,  a  great  portion  of  the  dose  received  in  the  basement 
from  the  field  near  the  model  came  from  radiation  penetrating  the  ground  and  con- 
crete.  This  is  pointed  out  in  Table  F-19  The  ground  penetration  was  determined 
using  the  same  procedure  as  that  used  for  the  larger  model.   However,  it  was 
necessary  to  take  ground  penetration  into  account  not  only  from  Area  I,  but  also 
from  a  portion  of  Area  II  (Fig.  9  ). 

The  raw  data  taken  on  the  1:12  model  begin  on  Page  77  of  Appendix  F  .   The 
normalized  data  are  given  in  Table  2. 
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Center 
of 

Model 


Area  I 

0.69  Sq.  Ft. 

l/8  Syrnme try- 


Portion  of 
Area  II 


Area  (in  ) 

Centroid 

Subarea  No. 

X 

-y 

1 

lU.07 

12. 45 

2.08 

2 

lU.07 

12. 45 

6.23 

3 

14.07 

12.67 

10.61 

k 

lU.07 

15.8U 

2.08 

5 

14.07 

15.84 

6.23 

6 

13.06 

15.28 

10.31 

7 

15.60 

19.07 

3.75 

8 

36.13 

22.87 

3.01 

9 

36.13 

31.31 

8.83 

10 

36.13 

18.30 

14.04 

11 

62.83 

29.83 

3.93 

12 

62.83 

27.80 

11.52 

13 

62.83 

23.87 

18.32 

Figure  9.   Area  I  of  the  1:12  model  field  and  the  portion  of  / 
from  which  ground  penetration  was  obt;< 
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SUMMARY 


The  normalized  dose  rates  in  the  last  column  of  Tables  1  and  2  give  the 
experimental  data  in  their  final  form.   The  standard  deviations  associated  with 
the  experimental  results  are,  for  the  most  part,  small.   The  largest  uncertain- 
ties occur  in  the  data  taken  in  the  basement  of  the  1:12  model  and  range  from 
7.37$  to  36.0$.   Elsewhere  the  standard  deviations  are  less  than  6,0h%   in  the 
1:12  model  and  do  not  at  any  position  exceed  7.56$  in  the  1:U  model. 

In  order  to  determine  the  protection  factors  for  the  various  positions  in 
the  two  models,  it  was  necessary  to  first  scale  up  the  model  data  to  correspond 
to  that  of  a  full-size  structure,  then  make  the  far-field  corrections  as  de- 
scribed  in  Appendix  D.   An  error  of  5%  was  assumed  to  be  made  in  performing 
these  operations.   The  protection  factors  found  in  Table  3  were  found  by 
dividing  the  scaled-up  corrected  numbers  into  the  free-field  dose  rate  at 
three  feet  (^485  +  11  r/hr  per  curie/ft  ).   As  shown  the  protection  factors  for 
the  two  models  agree  within  their  standard  deviations  except  in  the  corners  of 
the  basement. 

Also  found  in  Table  3  are  the  protection  factors  calculated  using  the 
"Engineering  Manual"  method  (Appendix  E)  and  those  determined  experimentally 
from  the  full-scale  block  house  during  summer  institutes.  Data  collected 
during  the  institutes  can  be  compared  with  the  average  model  results  in  five 
of  the  twelve  cases.   The  comparison  ranges  from  28$  in  the  center  to  29%  in 
the  corner  of  the  first  floor  and  73$  to  118$  for  the  same  positions  in  the 
basement.   Model  data  gave  the  higher  results  on  the  first  floor  with  the  op- 
posite being  true  in  the  basement. 

* 
Revisions  have  been  made  in  the  curves  used  in  reducing  the  model  data  sli 

the  above  results  were  obtained.   The  basem        lation.  to  be 

affected  the  rn< 
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The  engineering  method  was  used  to  calculate  the  protection  factors  for 
each  location  in  which  data  was  collected  in  the  models.   On  the  first  floor 
the  theoretical  calculations  were  conservative  "by  about  30$  when  compared  to 
the  model  data.   As  expected  the  two  compared  most  favorably  in  the  center 
position. 

The  engineering  method  has  consistently  failed  to  accurately  predict 

experimental  findings  in  basement  positions.   One  of  several  methods  suggested 

Ik 

for  improving  these  calculations  was  developed  by  R.  L.  French.    Unlike  the 

others,  however,  French's  method  requires  no  charts  additional  to  the  ones  in 
the  "Engineering  Manual"  and  is  still  comparable  in  accuracy  to  the  other 
methods.   The  engineering  method  basement  calculations  listed  in  Table  3  were 
obtained  with  the  aid  of  French's  method.   A  brief  outline  of  the  use  of  this 
method  is  given  in  Appendix  E.   After  the  adjustments,  theoretical  calculations 
were  within  25$  of  the  model  data  in  the  center  position  and  50%  in  the  corner 
when  a  portion  of  the  basement  walls  were  exposed.   The  comparison  dropped  to 
U3$  and  158$  for  the  same  positions  when  the  basement  was  fully  submerged. 

The  model  results  and  theoretical  calculations  show  that  the  safest  location 
in  the  structure  is  in  the  corners  of  the  basement.   The  results  also  show  that 
in  this  case  the  protection  factors  in  a  fully  submerged  basement  are  greater 
by  about  a  factor  of  five  than  those  for  a  basement  with  partially  exposed  walls. 

One  of  the  more  important  areas  for  improvement  lies  in  the  detection  system. 
It  is  expected  that  since  the  physical  size  of  some  of  the  dosimeters  used 
were  substantial  compared  to  the  dimensions  of  the  room,  the  measurements 
represent  an  average  dose  rate  over  a  volume,  rather  than  at  a  point.   New  de- 
tection systems,  such  as  thermoluminescent  dosimetry,  should  be  considered. 
Possibly  another  significant  factor  is  the  absorption  in  the  polyethylene  tubing 
between  the  model  and  source  of  radiation  emitted  in  a  direction  almost  parallel 
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to  the  ground,  several  effects  such  as  the  source  anisotropy,  source  energy 
degradation,  and  other  errors  introduced  by  the  approximation  of  a  point 
source  with  an  encapsulated  volumetric  source  probably  affected  the  results 
to  a  lesser  extent. 
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FACILITY  AND  EQUIPMENT 

The  experimental  portion  of  this  work  was  performed  at  the  Kansas  State 
University  Nuclear  Engineering  Shielding  Facility  (KSUNESF)  located  about  six 
miles  west  of  campus.   The  principal  facilities  are  as  indicated  in  Fig.  10 
and  include  a  full-scale  block  house,  foxhole,  5,000  square  foot  concrete  slab, 
1:12  and  l:k   scale  models  of  the  block  house.   The  equipment  available  includes 
a  pumped  source  fallout  simulation  system,  gamma-ray  projector,  sources  ranging 
up  to  about  80  curies,  fork  lift,  radiation  survey  meters,  dosimeters,  Gair   - 
Alarm,  and  other  instrumentation. 

The  concrete  block  house  from  which  the  models  were  designed  (Fig.  ll)  was 
constructed  of  lightweight  concrete  having  an  average  density  of  103  pounds  per 
cubic  foot.   The  roof  and  floor  were  made  of  pre-cast,  pre-stressed  concrete  slabs 
while  the  walls  were  concrete  blocks  stacked  without  mortar. 

Many  precautions  were  taken  to  assure  the  protection  of  all  personnel  from 
radiation  exposure.   The  experiments  were  performed  within  the  rules  and  regula- 
tions of  the  KSU  Radiation  Safety  Committee. 
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Figure  10.   Plan  of  shielding  facilities, 
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DATA  REDUCTION 

In  this  section,  a  resume  is  given  of  the  treatment  of  the  data  and 
associated  error  terms  from  the  readings  of  the  charger-readers  to  the  cal- 
culation of  the  protection  factors.   To  better  illustrate  the  techniques  in- 
volved in  reducing  the  data,  one  set  (Table  B-l)  will  be  taken  as  an  example 
and  all  operations  performed  on  it. 

Since  these  data  were  collected  with  2r  dosimeters  which  are  not  hermeti- 

o 
cally  sealed,  a  correction  was  applied  to  standardize  the  data  to  22  C  and 

760  mm  Hg.   The  correction  factor  (3  was  calculated  from  the  following 

expression: 

«    [273.0  j   0.555(T  -  32.0N1   760.0  ,R  ,x 

p  "        (295.0)(25.^0)  P  [B~1] 

where  T  is  the  temperature  in  degrees  Fahrenheit  and  P  is  the  atmospheric 
pressure  in  inches  of  mercury.   The  same  correction  was  made  on  the  data  ob- 
tained with  the  10  mr  dosimeters,  but  not  data  accumulated  with  200  mr  dosi- 
meters since  they  were  hermetically  sealed.   Table  B-2  gives  the  sample  data 
after  being  normalized  to  a  common  temperature  and  pressure.   Program  1  was 
used  to  perform  the  above  corrections. 

Since  the  exposure  doses  were  taken  in  terms  of  meter  readings  in  micro- 
amperes (\xa.) ,  it  was  necessary  to  relate  these  readings  to  the  true  dose  received. 
As  described  below,  a  calibration  line  was  developed  for  this  purpose  for  each 
dosimeter  with  the  aid  of  regression  analysis.   In  addition,  the  error  associated 
with  each  exposure  dose  was  determined  using  regression  analysis.   The  d 
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shown  in  Table  B-3  are  the  resultant  exposure  doses  and  error  limits  obtained 

from  the  calibration  lines.   Program  6  performed  this  function. 

o 
The  doses  were  next  normalized  to  r/hr  per  curie/ft  using  the  expression 


Dc  =  Cf  D>  (B"^ 


r/hr 

where  D   =   normalized  exposure  dose  m  h— — tttt? , 

c  curies/ft^' 

2 
A   =   area  of  contaminated  field  in  ft  , 

C   =   strength  of  the  source  in  curies, 
T   =   exposure  time  in  hours, 
and   D   =   dose  received  by  detector  in  roentgens. 
Table  B-H  shows  the  normalized  data  arranged  according  to  location  in  the  model. 

Since  these  data  were  taken  using  a  section  of  tubing  that  completely 
encircled  the  model,  the  normalized  doses  from  the  corners,  sides,  and  center 
were  simply  averaged  to  obtain  one  dose  for  each  of  the  three  detector  locations. 
In  other  instances,  the  simulated  fallout  field  covered  only  l/2  or  l/8  of  the 
perimeter  of  the  model.   A  different  averaging  procedure  was  followed  for  these 
cases.   First,  the  three,  four,  or  five  runs  that  were  made  were  averaged 
position  by  position  into  one  set  of  data.   Then,  for  the  l/2  symmetry  case  ,  the 
results  of  all  four  corner  measurements  were  added,  as  were  the  results  from  the 
side  locations,  and  the  sums  divided  by  two.   The  dose  received  in  the  center  of 
the  model  was  multiplied  by  two.   When  only  l/8  symmetry  was  used,  the  four  cor- 
ner values  and  side  dose  rates  were  added  and  multiplied  by  two.   The  dose  rate 
received  in  the  center  was,  of  course,  multiplied  by  eight.  Table  B-5 
summarizes  the  reduced  sample  set  of  data. 
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Each  of  the  other  sets  of  data  was  reduced  as  demonstrated  above.   The 
results  are  given  in  Table  1   for  the  l:k   model  and  Table  2  for  the  1:12 
model. 

The  normalized  doses  for  each  annulus  of  the  finite  contaminated  field 
were  added  for  each  model  configuration.   The  model  data  were  then  scaled  up 
to  the  full-size  structure  and  the  far-field  corrections  made  using  a  new 
method  developed  by  Kaplan  et.  al.  ,  which  is  described  in  Appendix  D.   The 
protection  factors  were  calculated  by  dividing  the  free-field  dose  rate  at 
three  feet  (U85  +  11  r/hr  per  curie/ft  )  by  the  experimental  doses  (Table  3). 

The  error  analysis  proceeded  as  follows:   Confidence  limits  for  each 
dose  were  obtained  with  the  dose  from  the  calibration  line.   When  the  ex- 
perimental measurements  were  normalized  to  r/hr  per  curie/ft  ,  the  error  terms 
were  also  normalized.   Since  the  source  strength  was  not  known  exactly,  the 
standard  deviation  associated  with  it  had  to  be  reckoned  with.   If  we  let  e-^ 
and  £2  represent  the  statistically  independent  errors  of  D  and  C  respective!;,', 
then  the  error  induced  in  Dc,  denoted  by  5,  as  a  result  of  the  errors  €j_  and 
e2  has  a  variance  equal  to 


If 


V(el)  +  Yi^A 
D2        C2 


(B-3) 


Sample  calculations  showed  that  the  errors  associated  with  the  exposure 
times  and  contaminated  areas  were  small  compared  to  the  errors  inherent  in 
the  source  strength  values  and  experimental  doses.  When  "n"  normalized 
doses  were  averaged,  the  variance  of  the  resultant  average  dose  was  found 

from 

V(Dl)  +  V(D2)  ♦  ...  4  V(Dn) 

V(D    )  «=  — •         ^l 

ave  P 

n 
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TABLE  B- 

1.   RAW  DATA 

(ua) 

Floor 

Elevated 

Area  II  (20". 
A=550  ft 

-160 

") 

1: 

:12 

Model 

First 

Floor 

s= 

=77. 

84 

Cur 

ies 

Run 

1 

Run 

2 

Run 

3 

Run 

;  4 

0)   <v 

6  e, 

t= 

=12.37 

min 

t=ll.98 

t= 

=12.25 

t=12.03 

co  3 
o  S 

co 
o 

T= 

=74°  F 

w 
0 

T=74 

co 
0 

T= 

=73 

co 
0 

T=76 

n 

P-, 

F-- 

=28.73 

in.  Hg 

Ph 

P-28.72 

Pn 

P=2 

Ph 

P=28.75 

1^7 

NW 

49.0 

NW 

49.5 

C 

47.0 

NW 

50.5 

134 

S 

42.5 

NE 

49.0 

N 

l*.5 

W 

4l.O 

137 

sw 

49.5 

SS 

U9.0 

W 

43.5 

S 

41.5 

135 

N 

^4-5.5 

C 

49.0 

SW 

50.5 

N 

45.5 

lUo 

W 

46.0 

E 

43.0 

NE 

53.0 

SE 

49.0 

151 

c 

49.0 

W 

44.0 

SE 

51.5 

C 

49.6 

150 

NE 

50.0 

SW 

50.5 

S 

45.0 

SW 

48.5 

145 

E 

kk.o 

N 

44.0 

NW 

50.0 

NE 

51.0 

lMl 

SE 

49.0 

s 

1*3.0 

E 

44.0 

E 

43.5 

TABLE  B-2.   DATA  NORMALIZED  TO  22°C  AND  76O  mm  Hg.  (ua) 


Floor  Elevated 
First  Floor 


Area  II  (20"-l60") 
A=550  ft 2 


1:12  Model 
S=77.84  Curies 


co 

CD 

Run 

1 

Run 

2 

Run 

3 

Run  4 

0 

Q 

CO 

CO 

co 

CO 

f£5 

O 

Pm 

t= 

=12. 

37 

min 

0 

Pn 

t= 

=11. 

98 

0 

Pn 

t= 

=12. 

25 

£  t=12. 

,03 

147 

NW 

51.3 

NW 

51.8 

c 

49.1 

NW 

53.0 

134 

S 

44.6 

NE 

51.3 

N 

46.5 

W 

43.0 

137 

SW 

51.8 

SE 

51.3 

W 

45.4 

S 

43.5 

135 

N 

47.6 

C 

51.3 

SW 

52.8 

N 

47.7 

i4o 

W 

48.1 

E 

45.0 

NE 

55.4 

SE 

51.4 

151 

C 

51.3 

W 

46.0 

SE 

53.8 

C 

52.0 

150 

NE 

52.3 

SW 

52.8 

S 

47.0 

SW 

50.9 

145 

E 

46.0 

N 

46.0 

NW 

52.2 

NE 

53.5 

144 

SE 

51.3 

S 

45.0 

E 

46.0 

E 

45.6 
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TABLE  B-3.   DATA  AS  INTERPRETED  FROM  CALIBRATION  LINES  (R) 


Floor  Elevated 

Area  II 

(20"-l60") 

1:12  Model 

First  Floor 

A=550  ft 

S=77.84  Curies 

Run  1 

Run  2 

Run  3 

Run  4 

t=12.37  min 

t=ll. 98 

t=12.25 

t=12.03 

1U7  NW  1.15  +  0.0738 

NW  1.16  +  0.07+1 

c  1.08  +  0.0726 

NW  1.19  +  0.0738 

13+  s  0.969+  0.0680 

NE  1.17  +  0.071U 

N  1.03  +  O.O689 

w  0.928+  0.0673 

137  sw  1.13  +  0.0793 

SE  1.11  +  0.0791 

w  0.956+  0.0767 

s  0.903+  0.0756 

135  N  1.04  +  0.0735 

c  1.13  +  0.0753 

sw  1.18  +  0.0761 

N  1.04  +  0.0735 

lHO  W  l.o2  +  0.0817 

E  0.9+2+  0.0807 

NE  1.21  +  0.0847 

SE  1.11  +  0.0830 

151  c  1.17  +  0.0778 

W  1.02  +  0.0753 

SE  1.24  +  0.0792 

c  1.19  +  0.0782 

150  NE  1.19  +  0.0722 

SW  1.20  +  0.0725 

s  i.o4  --  0.0695 

sw  1.15  +  0.071+ 

1I+5  E  0.993+  O.0767 

N  0.993+  O.0768 

NW  1.17  +  0.0795 

NE  1.20  +  0.0802 

lkk     SE  1.12  +  0.0830 

s  0.951+  0.0805 

E  0.978+  0.0808 

E  0.968+  0.0807 

TABLE  B-4.   DATA  NORMALIZED  TO  R/HR  PER  CURIE/FT2 


Floor 

Elevated 

Area  II  (20" 
A=550  ft 2 

-160") 

1:12  Model 

First 

Floor 

S=77.84  Curies 

Pos. 

Run  1 

Run  2 

Run  3 

Run  4 

S 

33.2  +  2.h3 

33.6  +  2.9+ 

36.0  +  2.52 

31.9  +  2.76 

E 

34.0  +  2.72 

33.+  +  2.9+- 

33.8  +  2.89 

3+.1  +  2.93 

N 

34.5  +  2.63 

35.1  +  2.82 

35.6  +  2.50 

36.7  +  2.70 

W 

35.1  +  2.90 

36.2  +  2.77 

33.1  +  2.7+ 

32.7  ■*  2.1+7 

SE 

38.5  +  2.96 

39.5  +  2.92 

1+3. 0  +  2.89 

39.1  +  3.0U 

SW 

38.7  +  2.8U 

1+2.5  +  2.71 

1+0.8  +  2.77 

ko.k  +   2.66 

NW 

39.2  +  2.66 

4l.l  +  2.77 

40.2  +  2.87 

42.0  +  2.78 

NE 

4o.6  +  2.62 

4l.2  +  2.67 

1+2.0  +  3.06 

1+2.1  +  2.96 

C 

40.2  +  2.80 

40.3  +  2.80 

37.5  +  2.61+ 

1+2.0  +  2.90 

TABLE  B-5.   AVERAGE  NORMALIZED  DATA  (R/HR  PER  CURIE/FT?)  FOR  THE  THREE 
DETECTOR  LOCATIONS 


Floor  Elevated 
First  Floor 


Location 


Side 

Corner 

Center 


Area  II  (20"-l60") 
A=550  ftP 


Norma]  d  zed  Data 


3k. k  +  0.681+ 
1+0.7  +  0.707 
1+0.0     5     I.39 


1:12  Model 
S- 77.81+  C 
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DOSIMETER  CALIBRATION 

Three  types  of  dosimeters  were  used  in  this  study.   The  Landsverk  2- 
roentgen  (L-8l)  dosimeters  were  used  whenever  possible  because  of  their  small 
physical  size.   As  small  as  they  were  (l.5"  long,  0.5"  dia.  ) ,  when  scaled  up  a 
factor  of  12,  they  would  represent  a  detector  1.5'  long  and  0.5'  diameter. 
Conversely,  the  Victoreen  10  milliroentgen  (mr)  dosimeters  were  used  only  when 
necessary  in  the  l:k   model  and  not  at  all  in  the  1:12  model  because  of  their 
size  (2.5"  long,  2"  dia.  ) .   A  Technical  Operations  charger-reader  was  used  for 
each  of  these  types  of  dosimeters.   The  third  and  most  frequently  utilized  type 
was  the  Bendix  200  mr  dosimeter.   A  Jordan  Electronics  portable  charger-reader 
was  used  with  these. 

For  each  of  the  2r  and  10  mr  dosimeters  a  calibration  line  was  developed  to 
determine  the  relationship  between  the  reading  in  microamperes  (\ia.)   on  the-  Techni- 
cal Operation  charger-reader  and  the  exposure  dose.   A  calibration  line  was  also 
developed  for  each  of  the  200  mr  dosimeters  since  the  reading  in  mr  did  not 
necessarily  correspond  to  the  true  exposure  dose. 

The  experimental  doses  were  obtained  by  exposing  the  detectors  to  a  known 
quantity  of  cobalt-60  radiation  for  varying  periods  of  time.   A  calibration  stand 
was  constructed  so  that  a  constant  source-to-detector  distance  could  be  maintained 
and  so  that  the  effect  of  scattering  from  the  ground  could  be  minimized. 

The  theoretical  dose  rates  were  calculated  using  the  following  equatio  : 

KSBBe  "^X 

D  -   ^4 ,  (C-l) 

k  it   X2 

in  which  D  =  exposure  dose  rate  in  mi lliroentgens  per  hou 

K  =  conversion  factor  I  (mr/hr)/(mc/cm2)~|  at  22  C,  760  mm  :  .  ;    '.  71  x  105  ^U^L2     , 


UO 


S  =  strength  of  the  cobalt-60  source  in  millicuries, 

18 
B  =  buildup  factor  due  to  air  scattering, 

Bp=  buildup  factor  due  to  ground  scattering, 

\x  =   total  gamma-ray  attenuation  coefficient  (cm  )  for  1.25  Mev.  gamma- 

-5   -l7 
rays  at  22  C,  760  mm  of  Hg;  =  6.79  x  10  J   cm   ,  and 

X  =  source  to  detector  distance  in  centimeters. 

As  an  example,  the  theoretical  dose  rate  for  the  200  mr  dosimeters  will  be 

calculated.   If  S.  is  the  strength  of  the  cobalt-60  source  in  millicuries  on  a 
l  to 

given  date,  t  is  the  time  in  years  from  the  date  corresponding  to  S.  to  the  date 
the  calibration  work  was  done,  and  t,  /?  is  the  half-life  of  cobalt-60  in  years, 
then  the  source  strength  on  the  date  of  the  experiment  is 

S  =  S±  exp  (-0.693  t/t1y2).  (C-  2) 

Substituting  in  the  appropriate  values  gives 


S  =  (265  mc)*  exp 
=  230  mc. 


-(0.693)(l.08)/(5.27v 


(C-3) 


*This  is  the  strength  of  the  source  on  August  10,  1965,  as  determined  from  a 
calibrated  Victoreen  Model  570-R  Meter.   The  other  sources  used  in  this  work 
were  calibrated  similarly. 


Ill 


The  buildup  factor  due  to  air  scattering  for  r  =  |-ix  =  6. 21  x  10   is 

Bx  .  i.o  *  (0.92)(6.21  x  10-3)  e0.0632(6.21  x  io-3)      ^_^ 
=  1.0057 

The  source  and  dosimeters  were  elevated  5'  2"  above  a  concrete  floor  and 
were  3'  0"  apart.   Using  these  factors  the  buildup  factor  due  to  ground  scatter- 
ing, Bp,  was  found  from  Clark  and  Batter  to  be  I.OO69. 

Finally, 

e~^x  =  e~r  =  0.99^ 
We  now  have  all  the  parameters  necessary  to  find  the  theoretical  dose  rate. 

D  (el)     =  (1.71  x  10^)(230)(l. 0057) (1.00686) (0.99*0  =  3?6  mr   (Q  _  ^ 
W  4(3.lUl6)(91.^)2  hr 

Similarly  for  the  10  mr  dosimeters  D  =  U9.U  ~   and  D  =  3.20  rp  for  the  2r 
dosimeters. 

The  theoretical  doses  were  plotted  on  the  abscissa  versus  the  experimental 
values  on  the  ordinate.   These  points  were  then  fitted  with  a  least  squares  line 
which,  for  the  200  mr  and  10  mr  dosimeters,  passed  through  the  origin.   For 
reasons  to  be  discussed  later,  two  lines,  one  of  which  passed  through  the  origin, 
were  necessary  to  describe  the  2r  calibration  lines.   Since  the  strength  of  each 
source  used  in  the  calibration  work  had  a  constant  percentage  error  associated 
with  it,  the  theoretical  doses  included  a  constant  error  tern.   Thus,  when  the 
least  squares  fit  of  the  data  or  regression  line  was  developed,  the  Berkson 
model   was  applied  to  take  into  account  the  constant  error  in  the  X  vail 
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Since  it  was  desired  to  determine  the  theoretical  dose  from  the  experi- 
mental reading,  a  regression  line  in  reverse  through  the  origin  (Eq.C-6)  was 
found  for  each  200  mr  and  10  mr  dosimeter. 

/      _,2  v   1/2 

1  ' 

where  y  is  the  mean  of  m  new  experimental  observations;  b  is  the  slope  of  the 
line;  $'  is  an  estimate  of  the  standard  deviation  including  the  error  associated 
with  the  calibration  source;  and  t  is  the  Student's  t-distribution.    For  a 

detailed  presentation  of  regression  analysis,  see  Statistical  Theory  and 

12 

Methodology  in  Science  and  Engineering  by  K.  A.  Brownlee.    The  constants  for 

the  regression  lines  for  the  200  mr  and  10  mr  dosimeters  are  tabulated  below. 
The  value  of  t  for  68  percent  confidence  (one  standard  deviation)  and  30  degrees 
of  freedom  is  1.01. 

The  10  mr  and  200  mr  dosimeters  were  essentially  linear  with  dose.   This 
was  not  the  case  for  the  2r  dosimeters.   The  calibration  curve  for  each  of  the 
2r  dosimeters  may  be  represented  by  two  regression  lines,  one  of  which  passes 
through  the  origin.   This  characteristic  is  a  consequence  of  the  dosimeter's 
construction.   Unlike  the  200  mr  and  10  mr  dosimeters  which  are  a  single  co- 
axial capacitor,  the  2r  detectors  have  a  center  electrode  made  in  two  sections 
which  are  not  in  contact  except  during  the  charging  and  reading  of  the  dosimeter. 
One  section  is  a  short  pin  centered  in  a  flexible  diaphram  at  one  end  of  the 
dosimeter.   When  the  dosimeter  is  being  charged  or  read,  this  end  of  the  dosi- 
meter is  pressed  into  the  charger-reader  receptacle  until  the  pin  makes  con- 
tact with  the  larger  central  electrode.   Both  electrodes  are  then  fully  charged. 
When  the  dosimeter  is  withdrawn,  the  pin  loses  contact  with  the  larger  electrode, 
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and  then  there  are  essentially  two  dosimeters  in  the  same  shell.   This  is  the 
reason  why  two  regression  lines,  rather  than  one,  are  necessary  for  the  . 
description  of  each  2r  calibration  curve. 

The  portion  of  the  curve  passing  through  zero  was  fitted  with  a  least 
squares  line  in  a  manner  similar  to  that  used  for  the  200  mr  and  10  mr  dosi- 
meters.  The  upper  portion  of  the  curve  was  represented  by  a  regression  line  in 
reverse  using  the  following  expression: 

*-  =  ♦¥  ^[(HV^fr^I72'  (c-7) 

where  k  is  the  number  of  observations  (x.,y.),   a  is  Zy./k,  and  the  other 
parameters  are  as  explained  earlier.   Again,  a  68  percent  confidence  level  was 
used  for  determining  t,  which  for  17-23  degrees  of  freedom  equals  1.02  (used 
with  Table  C- 3)  and  for  27-3^  degrees  of  freedom  equals  1.01  (used  with  Table  C-ty 

Tables  C- 3  and  C-U  list  the  constants  to  be  used  in  Equations  C-6  and  C-7  , 
respectively,  to  obtain  the  calibration  lines  for  each  dosimeter. 
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TABLE  C-l.      CONSTANTS  FOR  200  MR  DOSIMETER  REGRESSION  LINES 


Dosimeter  No. 

'  b 

s'2  -  (0, 

90 

1.071* 

2.358 

63 

1.029 

2.18U 

58 

1.0l*5 

1.116 

28 

0.979 

1.8M* 

1*8 

1.079 

0.721 

39 

1.059 

0.615 

62 

I.05I* 

1.107 

33 

1.016 

2.1U7 

70 

1.058 

0.653 

78 

I.089 

6.371* 

— 

I.0U5 

1.590 

86 

1.095 

1.969 

1*7 

I.058 

l.ooi* 

2lf (l/b2Ixn-2)   X106 


2.860 
3.121 
3.021* 
3.1*1*8 
2.838 

2.91*9 
2.972 
3.201 

2.953 
2.786 
3.021* 
2.758 
2.950 


TABLE  C-2. 

CONSTANTS 

FOR  10 

MR  DOSIMETER  REGRESSION  LINES 

Dosimeter  No. 

b 

s'2  -  (0.057)2        (: 

172 

10.  U8 

I.696 

168 

10. 1*8 

1.369 

165 

10.^3 

1.632 

176 

10.61 

2.338 

162 

10.80 

1.071* 

167 

10.6U 

2.558 

160 

10.56 

3.015 

166 

10.75 

1.738 

161 

10.66 

1.5^2 

(l/b2  £xj2)   X  lO^ 


1.557 
1.159 
1.168 
1.129 
1.1*68 
1.121 
l.ll+O 
1.099 
1.119 
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TABLE  C-3.       CONSTANTS   FOR  2R  DOSIMETER  REGRESSION   LINES    IN  REVERSE 
THROUGH  THE  ORIGIN 


Dos.    No. 


Useful  Range 


f2   -    (0.037)2  (l/b2  Ix,-2)   X  1QU 


lk7 

56.37 

y  <  29.9 

1.815 

131* 

52.17 

y(  31.2 

1.198 

137 

57.10 

y  <   29.5 

1.216 

135 

56.90 

y<  27.5 

I.U51 

l4o 

55.90 

y<   27.9 

2.864 

138 

56.94 

y<  26.7 

2.621 

151 

53.53 

y  <  29.6 

1.043 

150 

52.97 

y<  29.2 

1.126 

1^9 

52.65 

y<   30.  k 

1.355 

1U5 

57.79 

y<   26.2 

1.590 

1M 

56.74 

y<  28.8 

1.880 

1.394 
2.158 
1.358 
1.81k 
1.543 
1.913 
1.895 

2.147 
2.177 
1.413 
1.687 


TABLE  C-4.   CONSTANTS  FOR  2R  DOSIMETER  REGRESSION  LINES  IN  REVERSE 
(NOT  THROUGH  THE  ORIGIN ) 


Dos. 
Number 


Useful 
Range 


-n2 


b  s'2  -  (0.037)   1/k   b2  I(xj-x) 


lU? 
134 
137 
135 
140 
138 
151 
150 
149 

145 
144 


y  >29.9 
y>  31.2 
y>  29.5 
y>  27.5 
y  >  27. 9 
y>  26.7 
y>29.6 
y  >29.2 
y>  30.4 
y>  26.2 
y>  28.8 


1.115 

50.23 

34.91 

3.950 

0.0303 

7780 

1.099 

49.04 

35.35 

3.705 

0.0294 

8319 

1.115 

51.22 

36.41 

5.528 

0.0303 

8460 

1.099 

49.81 

36.02 

4.648 

0.0294 

8636 

1.099 

50.98 

38.26 

7.216 

0.0294 

9741 

1.099 

49.98 

36.74 

5.489 

0.0294 

8983 

1.115 

49.28 

35.15 

4.744 

0.0303 

7886 

1.131 

50.34 

36.22 

4.030 

0.0313 

7996 

1.094 

48.23 

34.69 

4.854 

0.0303 

7974 

1.099 

49.91 

36.55 

5.58.1 

0.0294 

8893 

1.099 

50.39 

36.30 

6.286 

0.0294 

8771 
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APPENDIX  D 
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NOMENCLATURE  USED  IN  THIS  APPENDIX 

DfX.w.h.r  — *■    r  ,  )   =  non-skyshine  component  of  the  dose  rate  in  a  structure 

v  '  '  '   n      n+1 

at  a  height  h,  from  an  annular  source  of  inner  radius 

r  and  outer  radius  r  , , . 
n  n+1 

h  =  detector  height  above  ground. 

Q  =  the  dose  rate  at  one  foot  from  a  one  curie  cobalt-60 

point  source. 
r  =  horizontal  separation  distance  between  point  source  and 

detector  in  the  structure. 
v.  =  inner  radius  of  finite  source  field  (effective  radius  of 

cleared  area  occupied  by  the  structure). 
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r  =  outer  radius  of  finite  source  field, 

o 

r  =  inner  radius  of  nth  source  annulus. 

n 

r  ,,  =  outer  radius  of  nth  source  annulus. 

n+1 

S(X,w,h,r  - — >-  r  ,)   =  skyshine  component  of  the  dose  rate  in  a  structure  at  a 

height  h,  from  an  annular  source  of  inner  radius  r  and 

'  n 

outer  radius  r 

n+1 

X  =  barrier  thickness  in  psf. 

a     (X,w,h,r)         =  skyshine  structure  attenuation  coefficient  at  a  height 

h  in  the  structure  from  a  ring  source  of  radius  r. 
A      •  =  scale  factor  used  for  obtaining  the  dimensions  of  a 

sca3e  model  from  the  corresponding  full-scale  structure, 
w  =  solid  angle  fraction  subtended  at  a  detector  by  a 

particular  barrier  in  a  structure. 
,j  =  attenuation  coefficient  for  cobalt-60  gamma  radiation 

In  air  at  22   C  and  760  mm  of  Hg. 
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SCALING  AND  FAR-FIELD  CORRECTIONS  TO  MODEL  DATA 

h 
A  new  and  significantly  better  method  has  been  developed  by  Kaplan,  et  al. , 

to  estimate  the  far-field  contribution  and  to  convert  scale-model  data  to  apply 

to  the  corresponding  full-scale  structure.   The  basis  of  this  new  method  is  the 

separate  treatment  of  the  skyshine  and  non-skyshine  components  of  the  dose  rate 

in  the  model.   This  procedure  is  necessary  because  air  density  is  not  scaled  in 

model  experiments.   Thus,  the  air  attenuation  in  the  model  experiment  is  different 

from  that  in  the  full-scale  case. 

In  order  to  use  this  new  method,  the  skyshine  and  non-skyshine  components 

of  each  of  the  measured  total  dose  rates  in  the  model  must  be  separated.   This  is 

done  by  first  obtaining  the  annular  source  skyshine  structure  attenuation  factor 

TO 

a        (X,w)  from  the  data  of  Burson  and  Summers   which  is  given  in  the  form  of  two 
s 

graphs  in  Kaplan's  work:   one  for  vertical  barriers,  and  one  for  horizontal 
barriers.   Therefore,  it  is  necessary  to  divide  a  structure  into  its  vertical  and 
horizontal  components,  find  the  skyshine  attenuation  factor  for  each  component, 
and  then  combine  the  factors  to  get  a      (X,w).   Having  this,  the  finite  field  sky- 
shine  dose  rate  may  be  found  from 

S   (X,w,hA,r.A— *r  A)  =  27rQk0»  (X,w)  e"^'A  -  e'^Vl.  (D -1  ) 

m   '  '  '  1     o        2  s  '  L  J 

This  result  is  subtracted  from  the  total  measured  dose  rate  to  obtain  the  non- 
skyshine  component  D   (X,w,hA,r.A — >-  r  A). 

The  experimental  data  may  be  scaled  annulus  by  annulus  or  over  the  entire 
finite  source  field.   Since  a     (X,w),  k,  and  k~  are  nearly  independent  of  r  over 
the  finite  field,  it  is  a  good  approximation  to  work  with  the  results  obtained 
from  the  entire  field,  rather  than  by  taking  each  annulus  separately  and  adding 
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the  results.  For  purposes  of  comparison,  several  sets  of  data  in  this  work 
were  treated  both  ways.   The  results  for  the  two  methods  were  almost  identical. 

The  expression  for  the  far-field  skyshine  contribution  in  the  full-scale 
case  in  terms  of  the  annular  shyshine  dose  rate  from  the  outer  annular  source 
in  the  scale  model  is 

SF.S.  (X>W>h'rN  —  ->  "  V^^N-l*  —  V>  >~ p^A   -WmA  "  (D  "2) 

e       -  e 

Therefore,  from  the  skyshine  dose  rates  in  the  model  for  the  annular 
sources  (calculated  from  Equation  D-l) ,  the  infinite-field  skyshine  dose  rate  in 
the  full-scale  structure  is 


e     -  e 


S„  -  (X,w,h,r1— *-  oo )  =  )   S  (X.w.hA.r  A  — »-r  ^nA) 
F.S. N  '  '  '  1     '        [_       m  '  '     '   n     n+1  '     -up  A    -up  ,A 

,  n       n+1 

n-1  e      -  e 

(  D  ^) 
+   S  (X.w.hA.r   a— ^r^) — r     ^  '^ 

mv  »  '      >   N-1      lTy   ~^PN-1     '^N^ 

e       -  e 

The  annular- source  non-skyshine  dose  rate  in  the  full-scale  structure  is 

f^1  5.  (up  )-E, (up  _) 

D    (X,W,h,ri  -*  -)  -  2   Dm(X,w,hA,rnA  -^r   A)   *   "   \     "+\ 

n=l  El(^nA)-El^Pn+lA) 


V^V 


+  D   (X,w,hA,rM  ,A — **r.T\)   =-? — = /\  P  / rr  ,   (D-1^ 

m   *  '  '  N-1     K      E-^np^Aj-E^up^)  » 
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where  the  first  term  on  the  right-hand  side  scales  the  model  data  up  to  full- 
scale,  and  the  second  term  makes  the  far-field  correction  to  the  data. 

The  method  described  thus  far  applies  only  to  above-ground  detector  locations. 
For  the  instances  when  the  detector  is  below  ground,  the  above  equations  still 
hold,  except  that  a  new  value  of  h,  the  detector  height  above  ground,  should  be 
used.   However,  little  error  will  be  incurred  when  r  is  substituted  for  p  in  the 
scaling  equations  (h=o)  if  the  building  height  above  ground  is  small  compared  to 
a  mean  free  path  in  air. 

To  illustrate  Kaplan's  method  a  specific  example  will  be  analyzed.   The 
case  will  be  considered  where  the  detector  is  in  the  center  of  the  basement  of  a 
one-story  rectangular  structure  shown  in  the  figure  below.   The  protection  factor 
will  be  found  using  the  data  from  the  1:U  model  of  the  structure. 


Figure  12. 
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Mass  thicknesses:   X  =  69  psf 

e 

Xf  =  55  psf 

X  =  55  psf 
r 

The  skyshine  attenuation  coefficient  a   (X,w)  will  be  examined  first.   The 

5 

ring-source  skyshine  attenuation  coefficient  through  the  roof  is  a        (X  +X  ,w') 

s  1     r     u  n 

(3) 
Next,   a        (X   ,w    )    /2     is  the  ring-source  skyshine  attenuation  coefficient   for 
s  e     3  V 

(2) 
the  first   floor  vertical  wall  at  the  center  of  the  floor  and  a        (X~,w   )      is 

s   v  f  u  H 

the  coefficient  for  the  first  floor  at  the  detector.   Therefore,  the  coefficient 

(3)  (2) 

through  the  floor  and  one  wall  is  the  product  a       (X  ,w  )-,,x  «   (X„,w  )TT   /  2. 

b    e  3  v  .  s  v  f  u  H 

The  total  structure  attenuation  coefficient  a   (X,w)  for  ring-source  skyshine 

s 


radiation  is 


aS3)(Xe'W^)VXffi2)(VWu)H 

SL)  fv    4.V      »•>  A       ),    5 e  ,    3    v s f u    H 


a  (X,w)  =  avx,(X_+X  ,v%     +     4  — - e      ^    '  s S-iL  (D-5) 

s     '  s         f     r'  u  H  2  7 

After  substituting  in  the  proper  values  and  using  the  graphs  in  reference  4, 
5rg(X,w)  =  </l}  (110,0.  24)R  +  2c/3)(69,0.20)v  x  ^2)  (55,0.  6o)R  .        (D-6) 

=  0.00112  +  2(0.032) (0.019) 

=  O.OO23I+ 
It  is  now  possible  to  calculate  the  annular  source  dose  rate  in  the  model 
using  Eq.  D-l  and  from  that  the  dose  rate  in  the  corresponding  full-scale 
structure  from  an  infinite  field  (Eq.  D-3). 

Sm(x,w,3A,ll.^?A  >  160/4)  -  2(3. Htl6) (13. 94) (0.U8) (0.00234)  x 

11.42         160 

(e"(WW  _  /(WW  } 

=  0.0073  r/hr  per  curie/ft 

Sm(x,w, 3/4,108/4, >-  160/4)  =  2(3.1^16) (13. 9>0(O.U8)(0.0023lO  x 

108  160 

( Q-(mm  _  e"(W8)(M  } 
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=  0.0025  r/hr  per  curie/ft 
Using  these  values  in  Eq.  D-3  the  infinite-field  sky shine  dose  rate  in  the  full- 
scale  structure  is 


11.1+2 


160 


^78                478 
SF^(X,w,3',H.U2  >co)   =    (0.0073)      6     lltk2  "     6 


160 


e'cm)w .  e~imiW) 


+  (0.0025) 


160 

WfB 


108  "  160 

(478) (4)  "(478)  (4) 


e 

2 


=0.096  r/hr  per  curie/ft' 
Next  the  non-skyshine  dose  rate  in  the  full-scale  structure  will  be 
found  using  Eq.  D-U.   The  non-skyshine  components  in  the  model  are  determined 
by  subtracting  the  skyshine  dose  rates  calculated  above  from  the  total  measured 
dose  rate.   Thus, 

D  (X,w, 3/4, 11. 42/4  >   160/4)  =  (1.53  -  0.007)  =  1.52  r/hr  per  curie/ft2 


and 


D  (X,w,3/4,  108/4  ^   160/4)  =  (O.O936  -  0.0025)  =  0.911  r/hr  per  curie/ft' 


m 


Using  this  in  Eq.  D-4  the  non-skyshine  dose  rate  in  the  full-scale  structure  is 


DF  g    (X,w,3',11.42' 


■*.«>)   =    (1.52) 


E. 


fll.421 


-     E. 


160 

L5rBJ 


'11.42  1       ~   r  160    -| 

(rnm)]     Ei  1(478)  (4)J 


E. 


+  (0.911)- 


160' 
LStEE 


] 


£5b 


it  r   100  1      „  r    lb0  n 

El[(W)(4)J        Ei  1(478)  (4)J 


=  1.60  r/hr  per  curie/ft 
The  total  dose  rate  to  the  fall-scale  building  from  an  infinite  field  is  1.70 
r/hr  per  curie/ft  .   Therefore,  the  protection  factor  is  - — =  or  286. 
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APPENDIX  E 
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THEORETICAL  PROTECTION  FACTORS 

The  protection  factor  is  defined  as  the  ratio  of  the  amount  of  radiation 
received  at  a  point  three  feet  above  an  infinite,  smooth,  uniformly  contaminated 

plane  to  the  dose  received  at  a  point  in  a  shelter.   The  methods  used  to  cal- 

c 

culate  the  protection  factors  are  described  in  Shelter  Design  and  Analysis, 

usually  referred  to  as  the  "Engineering  Manual".  This  textbook  was  intended 
to  aid  professional  engineers  and  architects  with  the  design  and  analysis  of 
structures  for  protection  against  radioactive  fallout. 

The  engineering  method  has  proven  reliable  for  calculations  involving 

17 
detectors  in  above-ground  locations,   but  it  is  deficient  in  the  analysis  of 

15 

basement  positions.     Several  authors  have  proposed  methods  to  adjust  the 

lk   15  16 
basement  calculations.   '  f  Each  of  these  methods  requires  additional 

charts  except  for  one  developed  by  R.  L.  French.   French's  technique  although 

just  an  extension  of  the  engineering  method  is  comparable  in  accuracy  to  the 

lk 

others.    The  example  below  illustrates  the  fundamentals  of  the  engineering 

method  and  includes  French's  extension. 

Given:   Structure  shown  in  Fig.  12. 

Detector  located  centrally  in  the  basement. 

Mass  thicknesses:   X  =  69  psf 

e 

Xf  =  55  psf 

X  =  55  psf 
r 

The  contribution  through  the  roof  and  vertical  walls  must  be  determined 
separately  then  combined  to  get  the  total.   In  the  extended  engineering  method 
the  ground  contribution  through  the  vertical  walls  is  composed  of  two  parts: 
That  which  comes  from  radiation  not  scattered  in  passing  through  the  floor, 
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C   ;  and  that  from  radiation  which  does  scatter  in  the  floor,  C   .   The  ex- 
pression  for  C   is 

cei  ■  vv^w^w-wl  L1- vvl 

♦[w-w]  W  E(e)}  t-vv] 


(E-l) 


where  the  notation  is  explained  below.   Had  the  non-barrier-scattered  fraction 

1  -  SW(X  )  been  omitted  from  the  above  equation,  the  formulation  would  have 
represented  the  original  engineering  method  expression  for  the  entire  ground 
contribution. 

The  contribution  from  the  floor-scattered  component  is 


C   =  B  (X  ,H)B'(X_)  <G  (w'^fl-S  (xj| 
g2    ee    of    auL    WeJ 


+  Gs(»;!SB(Xe)E(ejSs(VGb(»o)  . 

The  notation  in  the  above  equations  is  as  follows: 

B  (X  ,H)   Barrier  reduction  factor  for  exterior  wall  construction  as  a 
e  e 

function  of  the  wall  mass  thickness  and  detector  height. 
B1 (X  )     Barrier  reduction  factor  for  ground  contribution  for  floor 

immediately  over  the  detector  as  a  function  of  the  mass 

thickness  of  that  floor. 
G  (w  )     Accounts  for  the  skyshine  radiation  reaching  the  detector 

directly  and  also  that  reflected  from  the  ceiling  as  a 

function  of  the  upper  solid  angle  w  . 

G  (w^      Accounts  for  wall-scattered  radiation, 
s 

Sy(X  )     Fraction  of  emergent  radiation  scattered  in  wall  barrier. 

E(e)       Corrects  for  the  shape  of  the  building  for  v      utter 

W 
radiation;  a  fV     1  of  eccentricity  (e  -   -7—)  of  the  buil 

Li 

G,  (w  )     The  geometry  factor  giving  Lho  fraction  of  tl        ion 
b  u 
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scattered  in  the  floor  which  reaches  the  detector.   Assuming 
a  cosine  distribution  for  the  emergent  scattered  radiation 

G  (w  )  =  i_  (i  _  w  )2 
b  u  u 

w1 ' '       Solid  angle  fraction  subtended  by  the  walls  from  a  point  at 

the  center  of  the  first  floor  (not  shown  in  Fig.  12). 

Using  the  values  of  the  above  parameters  found  from  the  charts  in  the 

"Engineering  Manual",  the  equations  for  C   and  C   become 

SI      g2 

C   =  (0.19)  (0. 05^)  (O.09U-  0.079)(0.3U)  +  (0.U5  -  0.3k)  (0.66)  (l.klh)    (O.U) 
8l  '  L  J 

=  0.000M+  (E-  3) 

c   =  (o.i9)(o.05U^  Ro. 088)  (0.3U)  +  (0. 10(0. 66) (l.imol (0.6) (0.8*0 

=  0.00208  (E-  k) 

The  roof  contribution  C  (w",  X  )  is  a  function  of  the  solid  angle  sub- 

o  u'   o 

tended  by  the  roof  and  the  mass  thickness  above  the  detector.   This  contri- 
bution for  a  decontaminated  roof  was  found  to  be  0.000311. 

The  total  reduction  factor  R^  is  C  + C   +C   .   Thus, 

f  o       gl       g2 

R.  =  0.000311  +  0.0001^  +  0.00208  =  0.0028U,   and  P_  =  -~  -  353- 

I  I  JA_p 
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APPENDIX  F 
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TABLES  OF  DATA 

In  this  section  two  complete  sets  of  data  are  given:   l)  The  raw  data 
as  recorded  in  the  field  (Tables  F-l  through  F-3fr);  and,  2)  The  actual  doses 
and  their  standard  deviations  obtained  from  the  calibration  lines  (Tables  F-la 
through  F-3^a).   The  set  of  data  at  the  bottom  cf  each  page  is  the  data  at 
the  top  of  the  page  after  the  necessary  temperature  and  pressure  corrections 
have  been  made,  the  calibration  lines  applied,  and  the  data  rearranged  ac- 
cording to  location  in  the  model  during  exposure  (See  Fig.  6  ).   For  example 
in  Run  1  of  Table  F-l  the  reading  of  123  mr  from  dosimeter  number  58  be- 
comes 118  +  3.79  mr  in  Run  1  of  Table  F-la  in  the  NE  position  after  applica- 
tion of  the  calibration  line  for  dosimeter  number  58. 

The  symbols  t,  T,  and  P  in  the  data  represent  the  exposure  time  in 
minutes  (m) ,  outside  air  temperature  in  F,  and  pressure  in  inches  of  Hg, 
respectively.   It  should  be  remembered  that  the  temperature  and  pressure 
had  no  bearing  on  the  reduction  of  the  data  taken  with  200  mr  dosimeters 
although  this  information  is  supplied  on  the  tables. 

Each  dosimeter  is  numbered  according  to  its  type  as  shown  in  the  follow- 
ing table: 

TABLE  F.   DOSIMETER  FJMBERS  ACCORD IKG  TO  TYPE 

Dosimeter  Number  Type  Units  Read  From 

Charger-  Reader 

'  50-100  200  mr  Milliroentgen  (mr) 

100-160  2  r  Microamperes   (ua) 

l60-l80 10  mr Microamperes   (ua) 

The  units  of  the  data  given  in  Tables  F-l  through  F-3^-  may  be  found  using 

the  above  table.   The  units  of  Tables  F-la  through  F-3^-a  are  specified. 
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TABLE  F-la.   DATA  (MR)  AFTER  USE  OF  CALIBRATION  LINES 
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TABLE  F-2a.   DATA  (MR)  AFTER  USE  OF  CALIBRATION  LINES 
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TABLE  F-3a. 
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TABLE  F-k.      RAW  DATA 
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TABLE  F-l+a. 

DATA   (R)   AFTER  USE 

OF  CALIBRATION 
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TABLE  F-5a.      DATA   (IN  MR)   AFTER  USE  OF  CALIBRATION  LINES 
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TABLE  F-6a.      DATA    (MR)   AFTER  USE  OF  CALIBRATION  LUTES 
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TABLE  F-7a. 
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TABLE  F-8a.      DATA   (MR)   AFTER  USE   OF   CALIBRATION  LINES 
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TABLE  F-lOa.   DATA  (MR)  AFTER  USE  OF  CALIBRATION  LINES 


0 

•H 
-P 
•H 

W 
O 

Run  1 

Run  2 

Run  3 

Run  1+ 

FM 

t=l6.0m 

t=l6.0m 

t=l6.Qm 

t=l6.Qm 

88.9  +  2.81+ 

165  +  5.29 

1+5.9  +  1.86 


86.1  +  2.92 

156  +  1+.98 

1+1.7  +  1.50 


97.3  +  3.10 

160  +  5.20 

38.8  +  1.1+1 


91.8  +  2.93 

153  +  1+.98 

1+1+.1+  +  I.56 


69 


TABLE  F-ll.   RAW  DATA 
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TABLE  F-lla.   DATA  (MR)  AFTER  USE  OF  CALIBRATION  LINES 
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TABLE  F-12a.   DATA  (MR)  AFTER  USE  OF  CALIBRATION  LINES 
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TABLE  F-lUa.      DATA   (MR)   AFTER  USE  OF  CALIBRATION  LINES 
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Run  1                                Run  2                          Run  3  Run  h 

t=19.73m t=12.90m t=28.25m t^l9.05m 


S  53.6  +  2.0I+  37.2  +  1.8U  75.6  +  2.53  1+9.3  il.  81 

E  60.3  +  2.32  38.3  +  I.69  81+.  7  +  2.98  5I+.8  +  2.12 

N  52.6  +  1.90  35.0  +  1.80  73.1  +  2.59  50.2  +  2.11 

W  U6.3il.62  31.5  11.71+  63.2  +  2.1+2  1+3.1  1  1.79 

SE  80.9  +  2.88  51.111.71+  116  +  3.79  75.1  1  2.1+3 

SW  1+6.6  +  1.92  31.5  l  1.21+  65.1+  +  2.2U  39.1+  l  1.89 

NE  77.8  +  2.79  51.2  +  2.03  110  +  3.72  75.6  +  2.53 

NW  1+5.3  12.01  29.1+11.35  65.2  1  2.11+  38.9  11.87 

C  56.9  12.01  39.2  +  1.57  83.1+12.68  53.9+1.81 


lh 


TABLE  F-16.      RAW 

DATA 

Floor 

Elevated 

AREA  III    (27*    - 
A=1368  ft ^ 

ko') 

1 

:h  : 

Model 

Basement 

S= 

^76 

.17  Curies 

5-i 

-P 

Run 

1 

Run 

2 

Run 

3 

Run 

k 

•H   ,Q 

t=19.73m 

t=12. 90m 

_ 

t=28.25m 

m 

t=19.05m 

O     3 

0 

T=52 

w 
0 

T=31 

M 

0 

T=32 

w 
0 

T=22 

P    S 

Ph 

P=28.70 

Ph 

P=28.97 

CL, 

P=28.97 

Ph 

P-28.91 

166 

N 

3^ 

SE 

8.6 

E 

31 

NE 

13 

176 

WE 

12.5 

W 

29.5 

C 

89 

N 

3^.5 

162 

NW 

27 

NE 

9.5 

w 

65 

SE 

Ik 

161 

E 

21.5 

NW 

18 

NW 

39.5 

SW 

25 

172 

W 

kh 

SW 

35.5 

C 

62 

165 

SE 

13 

SW 

16.5 

N 

h9 

S 

32 

160 

S 

31.5 

c 

h3 

NE 

19 

E 

21 

168 

SW 

2U.5 

E 

15 

S 

ke 

NW 

27 

167 

N 

23.5 

SE 

20 

W 

J+5 

TABLE  F-l6a.      DATA    (MR)   AFTER  USE   OF  CALIBRATION  LINES 


0 

•H 
-P 

•H 
W 

Run  1 

Run  2 

Run  3 

Run  k 

& 

t=19.73m 

t=12.90m 

t=28.25m 

t=19.05 

s 

3.00  +  0.226 

i+.20  +  0.2^2 

2.88  +  0.192 

E 

2.03  +  0.157 

1.37  +  0.132 

2.76  +  0.187 

1.87  +  0.191 

N 

3.17  +  0.20U 

2.11  +  0.186 

I4J+9  +  0.262 

3.05  +  0.213 

W 

1+.21  +  0.251 

2.65  +  0.199 

5.75  +  0.31^ 

3.97  +  0.25I+ 

SE 

1.25   +  0.137 

0.763  +  0.13 

1.80  +  0.177 

1.22  +  0.115 

SW 

2.35  +  0.16J+ 

1.51  +  0.11+6 

3.2U  +  0.207 

2.20  +  0.162 

NE 

1.18  +  0.158 

0.839  +  0.106 

1.72  +  0.188 

1.13  +  0.136 

NW 

2.51  +  0.160 

1.61  +  O.lM 

3.5^  +  0.216 

2.U2  +  0.166 

C 

3.88  +  0.259 

8.02  +  0.356 

5.55  +  0.31*+ 

75 


TABLE  F-17.      RAW  DATA 


Floor 

Elevated 

AREA  IV 

(1+0'- 

60 

») 

1:1+ 

Model 

First 

Floor 

A=311+2 

ft2 

S=76 

.17 

Curies 

5-i 

Run 

1 

Run 

2 

Run 

3 

-P 

E    0) 

t= 

=67. 88m 

t= 

=58. 13m 

t= 

=  59-1+3m 

•H   ,Q 

w 

T= 

-.26 

• 

w 

T= 

=29 

W 

T= 

=31 

O    p 
Q  S 

o 

P= 

=29.05 

(J 
ft 

P= 

=29.05 

o 
ft 

P= 

=28.98 

62 

N 

75 

NW 

59 

W 

58 

90 

NE 

120 

C 

75 

SE 

103 

58 

E 

82 

S 

67 

NE 

98 

33 

SE 

111 

N 

61 

S 

61+ 

86 

S 

78 

SW 

61 

NW 

61 

63 

SW 

68 

E 

69 

N 

65 

— 

w 

70 

NE 

100 

SW 

61 

i+8 

NW 

70 

SE 

100 

E 

75 

70 

C 

82 

W 

59 

C 

72 

8 

•H 
-P 
•H 
W 


TABLE  F-17a.      DATA    (MR)   AFTER  USE   OF  CALIBRATION  LINES 


Run  1 
t=67.88m 


Run  2 
t=58.13m 


Run  3 
t=59.*+3m 


S 

E 

N 

W 

SE 

SW 

NE 

NW 

C 


71.3  +  2.51+ 
78.5  +  2.61 
71.1  +  2.1+0 

67.0  +  2.38 
109  +  3.69 

66.1  +  2.1+9 
112  +  3.75 

61+.  9  +  2.11+ 
77.5  +  2.50 


6l+.  1  +  2.21 
67.I  +  2.52 
60.0  +  2. 31+ 
55.8  +  I.87 
92.7  +  2.96 
55.7  +  2.11+ 

95.7  +  0.192 
56.O  +  1.98 

69.8  +  2.58 


63.O  +  2.1+2 

69.5  +  2.27 

63.2  +  2.1+2 

55.0  +  I.96 
95.8  +  3.29 
58.1+  +  2.16 
93.8  +  3.06 
55.7  +  2.11+ 

68.1  +  2.22 


76 


TABLE  F-18.   RAW  DATA 


Floor 

Elevated 

AREA 

IV 

(4o*- 

60') 

1:4  : 

Model 

Basement 

A=31)+2  ft2 

S=76 

.17 

Curl 

3S 

U 

-p 

Run 

1 

Run 

2 

Run 

3 

•H  ,0 

t= 

:67.88m 

• 

t=58.13m 

t-59.1+3m 

0  3 

w 
0 

T= 

:26 

0 

T=29 

m 
O 

T=31 

Q  S 

Ph 

P-29.05 

Ph 

F-29.05 

Ph 

P=28.98 

162 

SE 

19 

NE 

16.5 

C 

76 

172 

C 

90 

E 

25.5 

NE 

16 

166 

E 

29 

N 

4o 

SW 

30 

167 

NE 

19 

SW 

30.5 

E 

26 

165 

N 

45 

c 

71.5 

NW 

31.5 

168 

NW 

37 

SE 

15.5 

W 

51.5 

161 

W 

61 

S 

37 

SE 

16 

176 

SW 

35 

W 

53 

S 

39 

160 

S 

45 

NW 

33 

N 

4o 

§ 

•H 
-P 
•H 
W 

& 


TABLE  F-l8a.   DATA  (MR)  AFTER  USE  OF  CALIBRATION  LINES 


Run  1 
t=67.88m 


Run  2 
t=58.13m 


Run  3 
t=59.)+3m 


S 

E 

N 

W 

SE 

SW 

NE 

NW 

C 


4.01 
2.54 
4.06 

5.39 
1.66 

3.H 

1.68 

3.33 
8.08 


+  O.265 
+  0.179 
+  0.242 

+  0. 301 

+  0.128 
+  0.215 

+  0. 174 
+  0.204 

+  0.450 


3.29 
2.30 
3.52 

4.73 
l.4o 
2.71 
1.45 
2.96 
6.49 


+  0.205 
+  0.172 
+  0.218 
+  0.284 

+  0.133 
+  0.205 
+  0.122 
+  0.225 
+  0.360 


3.51 
2.33 
3.61 
4.69 
1.43 
2.66 
1.45 
2.88 
6.71 


+  0.231 
+  O.192 
+  0.249 
+  0.266 
+  0. 138 
+  0. 184 
+  0.146 
+  0.192 
+  0. 365 


77 


TABLE  F-19.   PAW 

DATA 

Ground 

Penetration 

AKEJ 

\   I  (-20") 

1: 

:12  iModel 

Floor 

Flush 

A=0. 

,69  ft2 

S= 

=0.227  Curies 

u 

CD 

-p 

Run 

1 

Run 

2 

Run 

3 

Run 

4 

CD  £-t 
S  CD 
•H  ,0 

t=42.22m 

t=42.22m 

t=>-!2.22m 

■ 

t-42.22ra 

o  3 

p 

T-36 

w 
o 

T-37 

0 

T-32 

w 
0 

T-35 

Q  S 

ft 

P=29.o4 

ft 

P-29.03 

Ph 

P=28.99 

Ph 

P=28.93 

, 

WE 

68 

H 

120 

WE 

70 

N 

118 

48 

N 

129 

NW 

64 

SW 

69.5 

ne 

65 

58 

WW 

68 

W 

84 

WW 

63 

w 

80 

39 

W 

79 

70 

sw 

72 

w 

82 

SW 

65 

90 

SW 

70.5 

63 

WE 

65 

w 

111 

nw 

54 

13^ 

SE 

40.5 

SE 

37 

SE 

36.5 

c 

15 

151 

E 

30 

C 

15 

S 

13.5 

SE 

31.5 

i4o 

E 

39 

E 

27.5 

S 

16.5 

135 

S 

16 

C 

15 

E 

27 

TABLE  F-19a.      DATA   (MR)   AFTER  USE  OF  CALIBRATION  LINES 


0 

•H 
•P 
•H 
CO 

Run  1 
t-42.22m 

Run  2 
t=42.22ra 

Run  3 
t=42.22m 

Run  4 
t=42.22m 

E 

N 

W 

SE 

SW 

WE 

NW 

C 


539  +  27.1 
120  +  3.80 
74.6  +  2.40 
813  +  65.6 
68.1  +  2.22 
65.I  +  2.33 
65.1  +  2.24 


271 

749 

115 

80.4 

721 

65.6 

63.2 

59.3 
270 


23.1 
78.7 
3.77 
2.67 
64.5 
2.47 
2.42 
1.98 
21.3 


241 

470 
108 

77.5 

698 

64.4 

67.O 

60.3 
252 


21.0 

35.1 
3.64 
2.50 
64.2 
2.12 

2.38 

2.11 
22.9 


284  + 

457  + 

113  + 
76.5  + 

576  + 
61.5  + 
60.3  + 
52.5  + 

277  + 


32.4 
26.7 

3.71 
2.56 
70.4 
2.03 
2.00 
2.16 
23.3 


78 


TABLE  F-20.   RAW  DATA 


Ground 

Penetration 

PORTION   OF  AREA 

11  (i 

20' 

•-3U») 

1:12  Model 

Floor 

Flush 

A-2, 

.06   ft2 

S=2.79  Curies 

u 

-P 

Run 

1 

Run 

2 

Run 

3 

Run 

1+ 

t= 

=1+9.  29m 

t=l+9.29m 

t= 

=1+9. 29m 

t=53.!+Ora 

o  3 

w 
0 

T. 

=50 

0 

T=56 

0 

T^ 

=58.5 

w 
0 

T=61+ 

Q  S 

p^ 

P- 

=28.60 

Ph 

P=28. 5h 

P-, 

P= 

=28.53 

a. 

P-28.1+8 

33 

E 

180 

W 

82 

C 

ikk 

86 

S 

132 

N 

60 

NE 

1+0 

C 

161 

1+8 

sw 

50 

E 

190 

N 

60 

S 

131+ 

70 

NW 

58 

S 

107 

E 

173 

NE 

1+0 

90 

N 

61 

NW 

57 

SW 

1+8 

— 

C 

150 

NE 

>+0 

NW 

58 

13J+ 

SE 

31 

11+0 

SE 

32 

SE 

31.5 

SE 

32 

E 

12.5 

58 

W 

89 

SW 

1+0 

S 

116 

63 

NE 

35 

C 

1^5 

W 

85 

TABLE  F-20a.   DATA  (MR)  AFTER  USE  OF  CALIBRATION  LINES 


0 

•H 

•P 

•H 

M 

O 

Run  1 

Run  2 

Run  3 

Run  1+ 

P-, 

t=l+9.  29m 

t=l+9.29m 

t=l+9.29ra 

t=53.1+Om 

S  121  +  3.97 

E  177  +  5.81+ 

N  56.8   +  2.26 

w  85.2  +  2.81 

SE  606  +  77.8 

SW  1+6.3  +  1.62 

NE  31+.0  +  1.78 

NW  5^.8  +  1.81+ 

c  ikk  +  1+.67 


101  +  3. 21 

176  +  5.66 
5I+.8  +  2.12 
80.7  +  2.88 

60I+  +  77.8 
38.3  +  1.55 
38.3  +  I.69 
53.0  +  2.17 

I'll  +  h.66 


111  +  3. 59 

161+  +  5.23 
55.6  +  I.87 
82.6  +  2.93 

621  +  77.9 
kk.7  +  1.98 
36.5  +  1.71 
55.5  +  2.09 

ll+2  +  1+.68 


12l+  +  3.9^ 

620  +  63.I+ 

37.8  +  1.39 

ll+7  +  1+.80 


79 


TABLE  F-21.      RAW  DATA 


Floor 

Flush 

AREA  I  (-20") 

1: 

;12  Model 

First 

Floor 

A=0, 

.69  ft^  . 

S= 

=0.227 

Cur 

ies 

<D 
-P 

Run 

1 

Run 

2 

Run 

3 

Run 

k 

0>  U 

S  0) 
i  <-*> 

t=i+2.22m 

t=U2.22m 

t=l*2.22m 

t=l+2.22m 

•H  A-1 

• 

to 

T=68 

w 

T=70 

w 

T-57 

w* 

T=5*+ 

O  p 

0 

a, 

P=28.92 

0 

P=28.89 

0 

Oh 

P=28.69 

0 

P-28.69 

135 

S 

20 

c 

22 

SE 

79 

E 

56 

1*5 

E 

*9 

s 

21 

C 

26 

SE 

70 

kS 

N 

170 

SW 

162 

w 

191 

NW 

128 

63 

W 

182 

NE 

105 

N 

155 

SW 

1^9 

1U7 

SE 

66.5 

E 

66.5 

S 

2U.5 

C 

22 

90 

SW 

160 

NW 

126 

NE 

107 

W 

182 

— 

NE 

109 

N 

165 

NW 

120 

N 

162 

58 

NW 

122 

W 

18k 

SW 

155 

NE 

98 

1U0 

C 

22.5 

SE 

69 

E 

57.5 

S 

22 

TABLE  F-21a. 

DATA  (MR)  AFTER  USE 

OF  CALIBRATION  LI 

NES 

% 

•H 
-P 
•H 
W 

0 

Run  1 

Run  2 

Run  3 

Run  h 

£ 

t=U2.22 

t=U2.22 

t=l+2.22 

t=*2.22 

s 

361  +  2k. 6 

375  +  25.6 

kkl   +  28.7 

397  +  33.7 

E 

1110  +  78.6 

i6ko  +85.8 

1290  +86.1 

1280  +78.2 

N 

158  +  5.03 

158  +  5.1U 

151  +  J+.97 

155  +  5.0U 

W 

177  +  5.83 

176  +  5.70 

177  +  5.70 

169  +  5.58 

SE 

1630  +  85.6 

1630  +  93.8 

19*0  +  96.5 

1660  +91.0 

SW 

l'+9  +  U.91 

150  +  U.79 

1U8  +  U.77 

1U5  +  J+.78 

NE 

10U  +  3.  *+5 

102  +  3. 1+7 

99.6  +  3.J+0 

93.8  +  3.06 

NW 

U7  +  3.76 

117  +  3.92 

115  +  3.77 

119  +  3.77 

C 

hlk   +  3U.0 

399  +  33.7 

U56  +  27.2 

393  +  27.9 

80 


TABLE  F-22.      EAW  DATA 


Floor 

Flush 

AREA  I    (-20") 

1: 

:12  Model 

Basement 

A=0, 

,69  ft2 

s= 

=0.227  Curies 

-P 

Run 

1 

Pain 

2 

Pain 

3 

Pain 

1+ 

•H  ,Q 

t=l+2.22m 

t=l+2.22m 

t=l+2.22m 

t=l+2.  22m 

o  3 

o 

T=68 

ra 
o 

T=70 

w 
0 

T=57 

w 
0 

T=5l+ 

«  s 

Ph 

P=28.92 

Pn 

P-28.89 

Ph 

P-28.69 

P-< 

P=28.69 

137 

S 

15.5 

C 

15.5 

SE 

1+1+.5 

E 

30 

151 

E 

29 

S 

15 

c 

17 

SE 

36 

28 

N 

125 

SW 

71+ 

w 

92 

NW 

59 

33 

W 

96 

NE 

70 

N 

125 

SW 

65 

13** 

SE 

35 

E 

30 

S 

20 

C 

86 

SW 

85 

NW 

69 

NE 

69 

w 

86 

70 

NE 

7U 

N 

133 

NW 

69 

N 

120 

39 

NW 

69 

W 

99 

SW 

83 

NE 

75 

11+1+ 

C 

17.5 

SE 

39 

E 

30 

S 

18 

§ 

•H 
-P 
•H 

w 

a 


TABLE  F-22a.      DATA   (MR)   AFTER  USE   OF  CALIBRATION  LINES 


Run  1  Run  2  Run  3  Run  1+ 

t=l+2.22m  t=l+2.22m  t=l+2.22m  t=l+2.22m 


Run  2 

Run  3 

t=l+2.22m 

t =1+2. 22m 

290  +  21.6 

389  +  25.3 

593  +  30.6 

5l+9  +  76.9 

126  +  3.99 

123  +  1+.10 

93.5  +  2.97 

9I+.O  +  3.22 

820  +  79.0 

9I+8  +  76.6 

75.6  +  2.71 

78.1+  +  2.52 

68.9  +  2.57 

63.0  +  2.32 

63.O  +  2.32 

65.2  +  2.11+ 

280  +  21.6 

322  +  22.2 

s  279  +  21.6 

E  560  +  70.2 

N  128  +  1+.23 

W  9k.  5  +  3. 26 

SE  729  +  61+.  6 

SW  77.7  +  2.71 

NE  70.0  +  2.27 

NW  65.2  +  2.13 

c  317  +  26.9 


320  +  27.0 

538  +  72.7 

113  +  3.60 

78.6  +  2.71+ 

7I+6  +  71.8 

61+.  0  +  2.1+1+ 

70.9  +  2.29 

60.3  +  2.32 


81 


TABLE  F- 

23. 

RAW  DATA 

Floor 

Flush 

AREA 

II    (2C 
ft* 

i"-l60 

") 

1: 

:12 

Model 

First 

Floor 

A= 

=  550 

S= 

=77. 

,8U 

Curies 

u 

Run 

1 

Run 

2 

Run 

3 

•H   ,Q 

• 

t=12.92 

t=ll.93 

t=12.08 

O    P 

03 
O 

T=7*+ 

w 
0 

T=78 

w 
0 

T=78 

Q  S 

ft 

R=28.7]+ 

ft 

P=28.72 

Oh 

P=28.71 

3.1*7 

s 

1+2 

SE 

1+7 

E 

1+2 

13U 

w 

»H 

E 

39 

N 

39.5 

137 

NW 

U9 

W 

1+0 

NW 

51 

135 

N 

1+7 

SW 

1+8.5 

W 

1+2.5 

3.1*0 

E 

1+6.5 

c 

1+7 

SW 

1+3.5 

151 

NE 

52.5 

N 

1+5 

SE 

1+8.5 

150 

SW 

52 

NE 

U8 

- 

NE 

50.5 

1*5 

SE 

50.5 

S 

1+1.5 

C 

1+1+ 

3M 

C 

1+7.5 

NW 

1+9 

S 

1+1 

0 


TABLE  F-23a.      DATA   (R)   AFTER  USE   OF  CALIBRATION  LINES 


Run  1 
t=12.92m 


Run  2 
t=11.93m 


Run  3 
t=12.08ra 


S 

E 

N 

W 

SE 

SW 

NE 

NW 

C 


0.931+  +  0.070 

1.01+  +  0.082 

1.08  +  0.071+ 

0.925   +  0.067 


1.79 
1.21+ 

1.27 

1.12 
1.08 


0.080 
0.073 
0.080 
0.080 
0.082 


0.930  +  0.076 
0.880  +  0.066 


1.06 
0.866 
1.09 
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TABLE  F-2l+a.   DATA  (MR)  AFTER  USE  OF  CALIBRATION  LINES 


S 

E 

N 

W 

SE 

SW 

NE 

NW 
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Run  1 
t=12.92m 


1+0.82  + 

1+1.71  + 

1+2.  5I+  + 

U3.O6  + 

1+7.1+6  + 
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TABLE  ?-25a.   DATA  (MR)  AFTER  USE  OF  CALIBRATION  LINES 
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TABLE  F-26a._     DATA   (MR)   AFTER  USE   OF  CALIBRATION  LINES 
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TABLE  F-27.   RAW  DATA 


Ground  Penetration 

AREA  I    (- 

■20") 

1:12  Model 

Floor 

Elevated 

A= 

=0.69  ft 

.2 

S=0. 

,227  Curies 

u 

(D 

-p 

Run 

1 

Run 

2 

Run 

3 

Run 

1+ 

Run 

5 

•H     rO 

t=70. 37m 

t= 

- 1+2.  22 

t=l+2.22 

t-- 

=1+2.22 

• 

t=l+2.22 

o  3 

o 

T=58 

0 

T= 

=56 

0 

T=56 

w 
0 

T= 

=56 

10 

0 

T-36 

Q   S 

Ph 

P=29.15 

ft 

F-28.68 

Ph 

P=28.67 

Pm 

P- 

=28.67 

Pm 

P=29.05 

33 

N 

118 

s 

175 

58 

W 

79 

C 

161 

N 

80 

C 

180 

1+8 

w 

50 

N 

80 

s 

181+ 

SW 

1+1+ 

86 

SW 

*+5 

70 

SW 

70 

N 

90 

W 

50 

C 

170 

S 

180 

39 

NW 

60-1/2 

NE 

81 

S 

163 

NW 

38 

w 

50 

63 

NE 

129 

SW 

1+0 

W 

50 

NE 

86 

90 

NW 

1+0 

NE 

83 

NE 

89 

NW 

39.5 

— 

C 

176 

NW 

38.5 

SW 

1+1 

N 

80 

151 

E 

55.5 

SE 

76 

E 

56.5 

13U 

SE 

73.5 

E 

57 

SE 

7l+ 

§ 

•H 
•P 
•H 
W 

a 


TABLE  F-27a.   DATA  (MR)  AFTER  USE  OF  CALIBRATION  LINES 
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TABLE  F-28.      RAW  DATA 
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TABLE  F-28a.      DATA   (MR)   AFTER  USE  OF  CALIBRATIOW  LIWES 
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TABLE  F-29a.      DATA   (MR)   AFTER  USE  OF  CALIBRATION  LINES 
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TABLE  F-30.      RAW  DATA 
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TABLE  F-30a. 

DATA  (MR)  AFTER  USE 

OF  CALIBRATION 

LINES 

0 

•H 
-P 
•H 
W 

Run  1 

Run  2 

Run  3 

Run  k 

O 

t=35.l8n 

t=U2.22 

t=U2.22 

t=U2.22 

S 

303  +  31.0 

373  +  2U.9 

397  +  33.7 

397  +  23.7 

E 

11U0  +  83.^ 

1220  +  79.8 

1260  +80.6 

1280  +78.I 

N 

llU  +  3.62 

133  +  ^.38 

132  +  1+.25 

135  +  J+.39 

W 

137  +  h.k5 

158  +  5.03 

153  +  5.03 

15U  +  If.  91 

SE 

1710  +  89.O 

i860  +  101 

1980  +99.6 

1880  +  101 

SW 

117  +  3.76 

139  +  Ml 

132  +  1+.19 

132  +  U.18 

NE 

81.2  +  2.86 

91.6  +  2.91 

96.  k  +   3.07 

93.8  +  3.06 

NW 

88.8  +  2.83 

108  +  3.50 

110  +  3.62 

103  +  3. ^9 

C 

301  +  21.9 

339  +  22.5 

316  +  23.9 

336  +  31.6 

89 


TABLE  F-31.      RAW  DATA 
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TABLE  F-31a. 

DATA   (R)   AFTER  USE 

OF  CALIBRATION  LINES 
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Run  1 

Run  2 

Run  3 

Run  1+ 
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t=12. 37m 

t=11.98m 

t=12.25m 
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s 

0.969  +  0.068 
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1.01+   +  O.O69 
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E 
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1.02  +  0.082 

1.02  +  0.075 

0.956  +  0.077 

0.928  +  0.067 

SE 

1.12  +  O.O83 

1.12  +  0.079 

1.21+   +  0.079 

1.11  +  0.083 

SW 

1.13  +  0.079 

1.20  +  0.072 

1.18   +  0.076 

1.15  +  0.071 

NE 

1.19  +  0.072 

1.16   +  0.071 

1.21  +  0.085 

1.20  +  0.080 

NW 

1.11+  +  0.071+ 

1.16  +  0.071+ 

1.16   +  0.080 

1.19  +  0.075 

C 

1.17   +  0.078 

1.11+   +  0.075 

1.08  +  0.073 

1.19   ?  0.078 

90 


TABLE  F-32.   RAW  DATA 


Floor 

Elevated 

AREA  : 

EI  (20"  - 

160  ") 

1: 

:12 

Model 

Basement 

A=550 

ft2 

s= 

=77. 

,8k   Curies 

-p 

Run 

1 

Run 

2 

Run 

3 

Run 

k 

a>  u 
g  a) 

t=12. 37m 

t= 

=11. 98m 

t=12.25m 

t=12.03m 

•H  rQ 

CO  g 

co 

T=7U 

co 

T= 

=7+ 

to 

T=73 

CO 

T=76 

O  p 
Pi  ^ 

O 

P=28.73 

O 

P= 

=28.72 

0 

P=28.71 

0 

P=28.75 

90 

¥ 

103 

S 

9+ 

SW 

109 

SW 

100 

63 

S 

95 

w 

92 

W 

100 

S 

97 

58 

SW 

100 

SE 

90 

c 

151 

w 

95 

1*8 

SE 

101 

E 

95 

NE 

112 

NW 

112 

39 

NE 

98 

SW 

96 

E 

101 

E 

98 

33 

N 

98 

NE 

9+ 

S 

96 

C 

1U5 

70 

NW 

102 

nw 

109 

N 

111 

N 

100 

— 

E 

100 

c 

1U0 

NW 

121 

SE 

102 

86 

C 

151 

N 

102 

SE 

106 

NE 

106 

TABLE  F-32a.   DATA  (MR)  AFTER  USE  OF  CALIBRATION  LINES 


>-< 

0 

•H 
-P 

•H 

CO 

Run  1 

Run  2 

Run  3 

Run  k 

s 

t=12.37m 

t=ll. 98m 

t=12.25m 

t=12.03m 

S  92.3  +  3.20  87.5  +  3.06  9^.5  +  3.26  9+.3  ±3.25 

E  95.7  +  3.19  88.1  +  2.82  95.^  +  3.03  92.6  +  2.9+ 

N  96.5  ±3.32  93.2  +  3.15  105  +  3.33  9+.5±3.01 

W  95.9  ±3.29  89.+  +  3.H  97.2  +  3.33  90. 9  +  2.98 

SE  93.6  +  2.99  86.1  +  2. 8U  96.9  +  3.26  97.6  +  3.25 

SW  95.7  +  3.12  90.7  +  2.89  101  +  3. +5  93.1  +  3.21 

NE  92.6  +  2. 3k  92.5  +  3.20  10^+3.30  96.9  +  3.26 

NW  96.  k   +  3.07  103  +  3.27  116  +  3.80  10H  +  3.30 

C  138  +  4.51  13*+  +  +.36  lkk   +  U.65  1U3  +  +.71 
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TABLE  P 

'-33.  : 

RAW  DATA 

Floor 

Elevated 

•  AREA  III  (160 

"-296") 

1:12  Model 

First 

Floor 

. A=1353 

ftd 

5=77.8*+  Curies 

Run 

1 

Run 

2 

Run 

3 

0>  U 

t=73.9m 

t=75.1m 

t=77.0m 

o  3 

w 
0 

T=8*+ 

w 
0 

T=6l 

0 

T=62 

P  S 

ft 

P=28.1+5 

£h 

P=28.55 

Ph 

P-28.53 

lU7 

W 

32 

E 

31.5 

NW 

35 

137 

c 

3^.5 

NW 

31 

NE 

32 

lH6 

E 

33.5 

C 

36.5 

S 

3h 

135 

sw 

33 

N 

31.5 

C 

37 

1U0 

NW 

35.5 

S 

33 

W 

32 

151 

S 

32 

SE 

3^ 

N 

32 

150 

SE 

3h 

SW 

35 

E 

30.5 

1^9 

N 

30 

w 

31 

SE 

35 

]M 

NE 

35.5 

NE 

35.5 

SW 

36 

§ 

•H 
-P 
•H 
W 
O 
P-, 


TABLE  F-33a.      DATA   (R)   AFTER  USE   OF  CALIBRATION  LINES 


Run  1 
t=73.9m 


Run  2 
t=75.3jn 


Run  3 
t=77.0ra 


S 

E 

N 

W 

SE 

SW 

NE 

NW 

C 


0.693 
0.718 
0.617 
0.662 
0.752 
0.689 
0.763 
0.765 
0.727 


+  0.071 
+  O.O69 
+  0.073 
+  0.067 
"  O.O65 
0.078 
0.078 
0.079 
0.07U 


0.676 
0.602 
0.612 
0.605 
0.706 
0.7^0 

0.715 
O.582 
0.758 


0.071 
0.066 
0.078 
0.073 
O.065 
0.073 
0.078 
0.073 
0.070 


0.688 
O.609 
O.652 
O.651 
0.720 
0.733 
0.613 
0.708 
O.762 


O.069 
0.072 
O.065 
0.071 
0.07*+ 
O.078 
0.073 
O.067 
0.079 
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TABLE 

F- 

■3h. 

RAW  DATA 

Floor 

Elevated 

AREA  III 

(l6o"-296") 

1: 

:12 

Model 

Basement 

A= 

=1353  ft 

■  2 

S= 

=77. 

81+ 

Curies 

u 

0) 

Run 

1 

Run 

2 

Run 

3 

-p 
0)    u 

•H   ,Q 

. 

t= 

=73.9 

t= 

=75.1 

t= 

=77.0 

CO    £ 

o  3 

o 

T= 

=81+ 

o 

T= 

=61 

W 

Q 

T= 

=62 

«  a 

fi 

P= 

=28.45 

Pi 

P= 

=28.  55 

PL. 

P= 

=28.53 

90 

SW 

18 

SE 

21 

N 

30 

63 

NE 

17 

C 

>+5 

NW 

18 

58 

S 

28 

N 

29 

E 

29 

1+8 

NW 

19 

SW 

19 

SE 

20 

39 

K 

30 

E 

29 

C 

>+9 

33 

C 

1+5 

WE 

18 

SW 

20 

70 

E 

25 

W 

29 

S 

30 

— 

W 

29 

s 

30 

NE 

19 

86 

SE 

18 

NW 

20 

W 

30 

PI 
o 

•H 
-P 
•H 

W 
O 


TABLE  F-3l+a.      DATA    (MR)   AFTER  USE   OF  CALIBRATION  LINES 


Run  1 
t=73.9m 


Run  2 
t=75.Im 


Run  3 
t=77.0m 


S 

E 

N 

W 

SE 

SW 

NE 

NW 

C 


26.8  +  1.31 

23.6  +  1.05 
28.3  +  l.lU 

27.7  +  1.1+8 
16.1+  +  1.39 

16.7  +  1.53 

16. 5  +  1.5*+ 

17.6  +  0.957 
1+1+.3  +  1.99 


28.7  +  1.50 
27.1+  +  1.12 
27.7  +  1.32 
27.1+  +  1.11+ 
19.  5  +  1.  56 

17.6  +  0.957 

17.7  +  1.55 
18.3  +  1.^1 
1+3.7  +  1.97 


28.1+  +  i.i6 

27.7  +  1.32 

27.9  +  1.67 

27.1+  +  1.51* 

18.5  +  0.973 

19.7  +  1.58 

18. 2  +  1.  3*+ 

17.5  +  1.55 

1+6.3'  +  1.59 
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APPEND  DC  G 


9h 


COMPUTER  PROGRAMS 

The  computer  programs  in  this  section  were  coded  in  Fortran  IV  for  use 
on  the  IBM  l'+Ol-l^lO  computer  located  on  the  campus  of  Kansas  State 
University.   Preceding  each  program  is  an  explanation  of  the  terms  used  in 
that  program.   The  reader  is  cautioned  that  a  variable  in  one  program  may 
or  may  not  be  designated  by  the  same  symbol  in  another.   Below  is  a  summary 
of  the  programs  used: 

1.  Dosimeter  Data  -  Normalizes  data  to  22°  C  and  760  mm  Hg. 

2.  Linear  Regression  I  -  Determines  a  calibration  line  through  the  data 
using  the  least  squares  method. 

3.  Linear  Regression  II  -  Determines  a  calibration  line  which  passes 
through  the  origin  using  the  least  squares  method  on  the  data. 

k.  Data  Interpretation  I  -  Uses  calibration  lines  to  calculate  actual  dose 
rates  from  experimental  data  from  10  mr  dosimeters;  also  gives  standard 
deviation  associated  with  each  dose  rate. 

5.  Data  Interpretation  II  -  Performs  same  function  as  number  k   but  treats 
data  from  200  mr  dosimeters. 

6.  Data  Interpretation  III  -  Interprets  data  from  2  r  dosimeters  which 
lie  to  the  right  of  the  intersection  of  the  two  lines  making  up  the 
calibration  curve  for  each  2  r  dosimeter;  also  calculates  standard 
deviations. 

7.  Data  Interpretation  IV  -  Same  as  number  6  but  treats  2  r  data  to  the 
left  of  the  intersection  of  the  two  lines. 

8.  Data  Normalization  -  Normalizes  data  to  r/hr  per  curie/ft  and  gives 
resultant  standard  deviations. 
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1.   Dosimeter  Data 

NOEXP  Number  of  experiments 

NPTEXP  Number  of  runs  per  experiment 

D  Measured  dose 

T  Temperature  in  F 

P  Pressure  in  inches  of  Hg 

TIME  Run  time  in  hours 

BETA  Factor  for  normalizing  data  to  2?  C  and  760  mm  Hg 


2?,  0 

203 

2n2 
203 

?  PC 


CORRECTIONS  FOR  TEMP.  AND  PRESS i,,  DETERMIna I  ION 
DlMENSI0ND(51>Vf(51)/P(51)»tlME(5l),AVGD6S(&l) 
FQRtfA i ( 1013) 
F0RMA1 (8F10.6) 

FORMAT (  ?UX.F10.6»4X.F6.2) 

forma i (  iox#6Kiu. 6#/J 


OF  AVERAGE  READI\S 


FORMAi  (  6Hl.3ETA  =  ,F10  .6) 
REAOU  >  2JH)  l^OEXP 
Do220J=li MOEXP 
REAU( 1,200 )NP fEXP 
IfTnPTEXP.EQ.O)    CAi.L    E  XiT 
READli,201)(D(I),l=l/LPTFXP) 

re :  a  d  ( i tZo'i  V  u 1 i ) » 1=1 »  nptexp  ) 

RF.Al)(x*20x)(p(T)#I  =  1,\3TEXP) 
RE  AD  1  1,201)1  I'-EU) 
XP  =  JPTFXP 
AVGDOS(J.)  =  0  .0 
D0218K=1, NP I EaP 


218 

?2_0 
221 


BET_As((273.  +  ,55S»J_T_<K)_-32ji  I  )*760  .  )  /  (  295  .  +  25  ,  4  01  #p  (  <)  > 

a'  R  I  T  E  I  3 ,  2  u  4  )  b  E  I  A 

D(K)=8t T A  *D  I  K  ) 

AVGDOSCJ)  -  A  v., D C b  (  j  )  *  (  D  (  K  )/XP) 

WR I T E I  3 , 2 1  3 ) ( D ( K ) i K s  J  , NPT£ XP ) 

WRI  rE(2#2u3>  (D(K)  »K»1#NPTEXP~J 

CONT r^uE 

D022l'I  =  li    0£XI 

WRITE(3,202)AVGD0S<l'#TlME(Xj 

no   ro  t>oo 

EnD 
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2.    Linear  Regression  I 

MODE      Number  of  sets  of  data  (one  set  for  each  dosimeter) 

Nl       Number  of  data  points  on  first  card 

N2        Sum  of  the  number  of  data  points  on  the  first  and  second  cards 

N3        Sum  of  the  number  of  data  points  on  the  first  three  cards 

NY       Total  number  of  (x,y)  values  in  each  set  of  data 

T        Value  of  Student ' s  t 

X        Calculated  dose  rate 

Y        Experimentally  determined  dose  rate 

The  remaining  terms  are  described  by  the  following  equation  -which 
was  used  on  page  k3   in  Appendix  C. 

XAVE  +  (Y  -  YINT ) /SLOPE  +  TOVB  *  SQRT  (SIGMA2  +  (0.03  *  Y)  **  2) 
*  SQET  (1.0  +  BUNK  +  (Y  -  YINT)  **  2/SIGMAD)  = 


y  +  y  1  a  +  tfs2  +   (Q.Q3y)2]  2 
x  ~T~~    -  b 


■1   1*      (y  -  a) 


2 


.1/2 


(m  +  k}  +  b?^(Xi-x)2 


:      LINEAR  REGRESSION!  ANALYSIS  FOR  STRAIGHT  LINE 
DIMENSION  X(35),Y(3d) 

10  0   FoRMAT(4I^,Flu.i>)  v 


101   FORMAT (  1QX,6F10.6> 

10  2   F  0 R M  A  T  (  /  ,  5 X  ,  F  1  0  ."  6  >  5  X  ,  F "  1  0  .  6  ,  5X7  F  1 0  ,  6  ,  5 X  i T 1  0  .  6  ,  5X ,  F  1  0  .  6  ,  5X  .  F  1. 5  .  9*  9 X  # 
IF  10, 6)  __  ____        . 

103  FORMAI  Cl3,"3F18.i5V 

104  FQRHATi  1UX»4HX  =  , Fll . 6 , 4H+0R- , Fll , 6     ) 
1 0  5  T 0 R M A  T (8  F l  DTZT " 

106  F0RMAT(2I3) 

107  FORMAT  ('  1UX,4F10  .6) 
23r   FO^MAl (5Fi0.6,Fi5.9) 

99   FORMA'T"('83hX=  XAVE  ♦  (Y  •  YtKTT) /'SLOPE  tOR- ~ TO VB*SPRIME I  1 /M  ♦  BONK  ♦ 

1  ( Y  -  Y J N J  >  *  *  2 /  SIGMAD>**i5) 
9fi   FORMATi//,5Xf«HXAVF*llX,4HYlN'T,llX,5MSL.  OPE#10X#4HTOV8#11X»4HBUMK#1" 
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11 XV  6  HS I G  M  A  D~«  1 4X#  &H5  £  6H  4 1) 
204       READU,10U)MODE 

IF  (MODE.Gl.  0)    GO    TO      203 
SUMX=u . 

S  U  M  Y  =  U  . 

SUMXYsO  . 

S  U  M  Y  2  -  0  , 

S  u  M  X  2  -  0  . 

R£~AD0.4  1 :-'  'j  JNl",  N2,N3,  NY*  T 

M  11  =  N  1  ♦  1 

\j22  =  \2  +  l 

N33  =  N'3  +  1 
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READ<_1'*01J 
R  £  A  D  < 1 » 1 0 1) 

R£AD(1j101? 

r  e  a  b  <  i  #  i  n  l ) 

REA0^1»101J 

READflViOi) 

R  E  *  D  (  1 >  1  0  i  ) 

READTi"»ioi) 

DO_200lfl«NV 

SUMX  =  SUMX"i*X(i) 

S  y  M  Y  =  S  U  M  Y* Y  ^1 J 

SUMXYsSUMXY*  X(i")+Y(I) 


IX(I),T=1»N1) 
I  X  (  J  )  ,  J  =  k  1 1 ,  m  2  ) 

(X(K),K=N22*N3) 
l  xn. )  ,UsK-33,NY) 
|  Y(  T)  »  l  =  1»M) 
i  Y<  J)  ,  J  =  Mll,N2>" 
CY(K)  >K  =  ,\|22>N3) 
(Y(L)il=N33»NY) 


2on 


203 


SUM 
SUM 

FY  = 

Su'O 

XAV 

YJM 


S  U  M  '  ?  +  Y  (  I  )  *  Y  <  I  ) 
=  SU  1X?*X(I)*X(I) 


Y2 
X? 
NY 
K« 

PE 
E  = 
r= 

[GMAD=SL0PE**2*(SUMX2-SUMX**2/FY) 

MA 
MA 
B  = 


1 .  /  F  r 

=  (S'J;:xV-5UmX*SUmy/f  Y  j  /  (  SUMX2-  SUmX*  *2/F T  j 
:  S  U  M  X  /  F  Y 
iSUM"Y7f  f 


SIG 

Sir, 

TOV 
WRi 


2=(Somy2~SUMy**2/FY-SL0PE**2*<SUMX2 
1=SQRT (S1GMA2) 

1/SLuPE 


SU^X**2/EY)  ) / (FY-?,  ) 


(3,99) 


WRITE (3. 98) 

W  R I T  E  f3 » 1 0  2  TX  A7  E ,  Y I N  TT^L  0  P" E  #  TOVI  #  B  U  N  K  #  S IG  MTD  #  S 1 G  M  A 1 

WRITE (2* 230 )XAVE i YpNT# SLOPE* TOVB* BUNK, SIgMAP 

GO  TO  204 

IF (MODE  .EU.  J  )CALL  EXIT 

E\D 
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3.    Linear  Regression  II 

The  following  equation  defines  the  terms  that  were  either  changed  or 
not  used  in  the  preceding  program: 

Y/SLOPE  +  SIGMA  *  SQRT  (l.O  +  SIGMA1  *  Y  **  2)  = 

1/2 


y  +  ts 
b     b 


2     n 
i     y 
+  b^B^iT 


m 


NOTE:      The   constant  percentage   error   in  the  x  values   due  to  the   standard 
deviation  associated  with  the   calibration   source   strength  was   included 
in  the  error  term  above  by  hand. 


10P 
101 
102 

103 
10.4 
105 

106 
204 


20  0 


203 


LINEAR    REGRESSION    ANAlYsIS    FOR    STRAIGHT    LIVE    THRu    ORIGIN 

DIM£NSIQN    X(4Q)»Y(4Q),TIM£(40)  

FqrmaT (    I3,F1Q.5) 

F  0  R  M  A  T (     1  u  X  ,  3  F 1 0  .  6  > 

FORMAT!    10X,6hX    =    Y>,Fll.6"#4H+OR-#F10"l8/8H(l/M*Y2*F18t15»lH)) 

FORMAT (13, 3F18.15J 

FORMA  T(  10"X'V"4HX  =  iF11.6*4H*0R-iFll,6     ) 

FQRMA.TJ  oFlG , 6) 

F  0  R  M  A  JJ  1  •;  X  ,  F 1 1  .  6  ,  F  U  .  3 ,  F 1  8  , 1 5  > 
READ(l»i0  0VM6DE 
IF(MODE.GT.O)  GO  TO   203 
suMXY=n . 

SUMY?=C,  

SUMX2T6, 

READU,10  0)NY,T 

XNj  =  NY 

RE  AD (1,101)   (X(I)  ,  E  =  1,NY) 

READ(l,10b)  ""'"(  YtI),T  =  l,NY) 

Dp  2  0  U  1=1, MY  

SUMXY  =  SUM~xY+  X(I)*Y(I) 


SUMY2 
SU^X2 
SlOPF. 
RAT  =  < 
SIGMA 
SIGMA 
WRITE 
WRITE 


=SUMY2+Y(I)*Y(I) 

=SUMX?+X( I)+X(I> 
=SUMXY/SUMX2 


SUMY2»<  (SUMXY*SU'1XY)/SUMX2)  )/  (XN  »1.  ) 

?"<"sqkt"(Rat  >*♦  r )  /slope 
i  =  i. /(Slope* Sl op e*sumX2) 

li,lU?)SLOrE,SiO-'A,SlGMAl 

(2,l'i6  )Sl0kF,RAT,SigMA3l_ 
GO  TO  2Q4 
IF(M0PE.£Q.1>CAIL  EXIT 

END 
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h.        Data  Interpretation  I 

NOEXP  Total  number  of  sets  of  data 

NPTS  The  number  of  repetitions  of  an  experiment 

L  Dosimeter  number 

SET  The  number  of  the   set   of  data 

X(L)  A  calculated  dose  rate  for  a  particular  dosimeter 

Y(L)  An  experimentally  determined  dose  rate  for  a  given  dosimeter 

SIGMA  The   standard  deviation  associated  with  each  data  point 


1.P8 
HO 
U2 


[NT 

TOR 

ro* 

113) 

moo   «ea 
IF< 

00^ 
REA 
REA 


ERPFETATIOiM    OF    DATA    FUQ"    10       MR    D0SIME1FRS 

ENSjON    Y(?^>0) 

M  A J  (  2 1 3  ) 

MAT( 1Qa,^F10. 6) 

MAT  ("/f  10Xri"HX#I3#3H    =    ,F10.6,5X,/H+ 


OR 


#F10  ,6,bX,9HDATA  SE1  » 


D( 1, 108 > NOEXP 
NOEXP. LB. 0  )CALtExlT 
Q0J=1. ^TEXH 
D(1/108)NPTS 
DTa.,io6>L>SEl 


R"EAD<1»U  U5  (V(T)/I  =  1»kPTS) 

D03QOI  =  1,JvPTS 

TF  ti.EOY!  /?)G0TU1Q 

EG.  168 )GOTO20 

EC  .  16?  )  GOT  1)3  0 

E  G . 17  6  >  G  0  T  0  4  0 

EC  .  1  62>(iO~fu5Q 

EO.167>GOTO60 

E  &  n  "6  0  )  CTO  T  0  7  0 

■      .  1  o  ft  )  G  U  T  0  8  0 

t  C  .  1  o  1  )  G  u  T  u  9  0 


in 


IF(L 
IF(L 

if  a 
iF'a 

[F(L 

ifTl 

IFJL 

if  a 

X172 


v  (  I  i/i  1  .3.54285 


SIGMA  s(  l.Q  1/H.534265  )*SQRT(1«  69598*  ( ,05»Y(I) )  »♦?  )  »SQRT(1  .♦  t  00  0  013 
■   -  r046'2964*Y  <T)*V2j 
dRITE(3, 112  )L,X172, SIGMA, SET 
G  0 1 0  3  0  0 
20    xi>8=    MI)/    11^316036 

S I GM A  s  ("i .  Vl  / 1 1 V3 1 6 0 36 )  ♦  SOR T*(  i 
1   i    &69821*Y(I)**2) 
WRlTE(3,112)l»~Xl6    ■     '  }Ma#SE1 
rOvSOO 


368/4*(,ni>*Y(l))**?)*Si:RT(i.*,uO0Ull 
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3   X165  =  Y<I)/il.269/12 

SIGMA =( 1.0 1/11. 269712 >*SQRT(1.63l69»( t 05* Y( I) )**2> ♦SQRT<1 .», 000U11 


WRIT  E(3V112TL'X165#"SIGHA»  SET 

G0T03QQ 
40  X 1 7  6  =  V(l)/ll. 402611 

SIGMA=(1. 01/11. 4626li)*SQRT(2,33798*(i05*Y(I>>**2>*5GRT(l. +.000011 
1292268393*Y(I)**2) 

WRITE (3, 112 )L»X176, SIGMA* SET 

G0fo30d 
50  XI 6 2=  Y(I)/11.6/20  0  9 

S1GMA=(1.  i;i/ 11.672009  )*SQRTC  1,074  ij9+(  ,  05*  Y(I)  )**?)*SORT(l.-t-,  000014 

1682449470*Y(I)**2) 

w  R  1 1  E  (  3  ,  1  i  2  )  L  ,  X 1 6  2  ,  5 I G  *  A  ,  S '  E  f 
GOTO300 


Ail  Xl6/=  Y(I) /U  .5U2398 

SIGMA  =  ( 1 .  Ul/ll  ."&02398 )  *SORT  ( 2 .  55768* (  ,  05  *y  ( I  ) )  **2  >  *SQRT  ( 1 .  * .  000  Oil 
1214283768*Y^IJ**2) 
Wr ITE  (  3 7  112" ) L »  X167 *  SIGMA , SET 

GOTU300 

70  X160"=  Yd  )/ll,4~0W69 

5 IGMA= ( 1. 01/11. 408569 )*gQRT( 3.0 1494* ( .  Q5*Y(I) > »*2 ) »SQRT ( 1 « +  . 00  0  Oil 
1399504586* Y< I  > **2> 
WRITE (3,1 12 )L»X16  0, SIQMA »SET 

GOTO3O0 
RO  X166=  Y(I)/11. 617152 

SIGHA  =  (l.1ji/ll.t>l/152)*sf^T(1.7  3  806*(.05*Y(I))**2)*SQRT(l.  +  .000U10 


1993828 15 0*Y( I) **2) 

WRITE(3>112)L»X166#SIGMA>SET 

GOTO3U0 
go  X161?  Y(I)/11. 516460 

SIGMA=< 1. Ul/ll. 516460)+S OR T(l . 5  42  4  5* (  ,  ~05*Y(1  )  )**?)*SOHT(3  .  +  .0  0  0011 
1166913951*  Y(I_)_**2) 

WRITE  ("3#  112)  L»XX61#SIGMA/SET 

GOT  03 0  0 

300  continue! 

40  0  CONTINUE  

GGTQ1U00 
END 


101 


5.   Data  Interpretation  II 

The  meaning  of  the  terms  in  this  program  are  the  same  as  in  program  h, 
DIMENSION  Yllu) 


1  „ p  rn^M^fPT^ 
1 C  9  FOR  M  A  l  ( b>  r  1  0  .  6  ) 

112  FORMAT  (7 7l(fX»~iHX«l3*3H  =  »F£0 .6*5X» 7H*  OR  -  ,  Fl  0  .  6  » >5  X  ,  9HOATA  St  I  » 
113) 

INTERPRETATION  OF  DATA  FRO'^  200  Mfl  DOS ItfFTFNS 
10  0  0  R£AD<l'l08)f<o£XP 

IF ( NOEXP.Efl. 0  >CALIEXIT 

D04UOJ=1,  OEXP 

R  E  A  D '( 1  /l  0  8  >  N  P  T  S 

READ(I#108)L»SE1 

RLaD<  1*  109)  (  Y(  I;  ilsiVNPf'S'j 

D0300Isl»NPTS 


[F(U_.  EC.  90)  GOT  O10 
I F  ( U  •  E  6 .  6  3  >  G  0  T  0  2  0 
IF(L.EQ,58)G01030 

I F  (  L  •  E  Q  .  2  6  )  G  0  1  0  4  0 
TF ( L, K 0.4m) gq, Qb 0 
I FTL7E07  3  9TG  07060 
IF(U.£O.62)G0 • O/O 
IF (U.EO.33)GO"O80 
IFCU.EQ, 7  j ) G G  J  0 V 0 
IF  f l~.EC . 76>G0 >  0110 
IF(U«£Q»00)GOTO130 
IF  f  L-  Eti.86)GO'iO150 


[F(|.tE0t  47)601  0160 
1 0  X 9  0  =  Y  <  I )  / 1 . 1 6 1 23 7 

SIG"MAsri."0l/iTi61237T*SQRT(2.35794*C.03*Y(I)  )  **2)  *SQRT  (1,*.  OOOOOZfl 
1 6  0  3  9  8  0  4  2  *  r  ( 1U  *  2 ) 
WRlf£(3VlY2>l/X90#  SIGMA*  SET 

Goro30  0  

20  X63VY('i)/l".  111735 

SiGMA=<l.'-l/l.lll/35>«SQRT<2.1837;>*(.03*Y(T>)**2)»StiR'<l.  +  .OOOOU31 
l20794!>6l*Y(I)**2) 

HRITE(3,112)L*X63»SIGMA<SET 
G"oTO3"00 

30  X5B*Y(I)/i, 129400 

SIGMA* (l.Ul/1.1294QQ)»SQRT(l. 11626* (  ,03*Y<D  )  »»2)  ♦SQRT(1 1 ».  Q00003U 

I  2  3  9  J  4  0  4  9  *  Y  (  I  )  *  *  2  ) 

MRITE(3»112)l.#X!>8j»SIGMAjSET 
QOT03 
40     X28sY(I>/l,057666 

SigMA  a  (1\  u'l/d  ,  Oi>7666)*~3QRTUY8439T*l,oy*Y<:J~)  iTt*2T*S"QR"l  (l.*,orooj34 
148025492*Y(I>*#2) 
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WRI  r£(3,  112)L»X28'SXGMA7SET 

GOTQ300  

br,    X  4  8  =  yY'j  >  / 1 . 1 6"F8  0  4 

S  IGMAsd.jjl/ 1.16580  4)  *SQRT(     .72Q77+(.Q3»Y(I) )  **2 )  »SQRT  (1 .  + 1  Q0Q0028 
138 0  28652*^  (T>**2) 


WRI TE 1 3 , 112 ) L  * X48 »  SIGM A , SET 

GOT  03 0  0 

6  n  x  3  9  =  Y  ("I  >  /~1 ,  14  3"7  0  6 

SI6MA=(1.U111.1«37Q6)*SQRT  (  .  63  45  0  +  (  ,  Q3_*  Y  ( I  > )  *  +  2 )  *SQR1  (1,  +  ,  00000 29 
1 4 8 7 56 3 1 0"* Y~(  I T* V2 ) 
WRITEC3»li_2)L.X39»SIGMA,StT 
GQTQ300 

7  0  X62  =  YJI>/1. 139273 

SIGMA=(1.  G  1/1.  13  92  73  )+SQRT(i.  10/0  /+(  .  03*Y(i)  >**2>'*SQR  \  f  1.  +  .  0000029 
1  7174863  O^Of  J  11**2) 
WRI TET'3j 112  )L  iX62/BTiGMAVS"Ef 
G  0  T  0  3  0  0 


8  0  X33-V ( I )/i. 097697 

SlG"MA=iTl."uT/l.  0  V  769  7  >*SQRT(  2.' 1466  (]  +  ('■;  0  3+Y<I>  )**2)*SQPT(l.+.O0OGG32 
101124866»Y(I)**2) 
WRI  rE( 3/112  )LiX3"3, SIGMA, SET 
GOTO30  0  _  

9  0  X70=t  (I  )/1.14"29o7 

SIGHA=(1.U 1/1. 142967 )»SQRT<  . 6  534  5+ ( . Q3+Y(I) ) **2) »SQRT ( 1 .+ . 0000  029 
1525698'26+Y"d5*"*2 ) 

WRITE(3,ll2?L,X7Q,SIGMA,S£T   _____ 

GOTO2T00 
110  X78-VH  )/i.  176685 

SIGMA=(1  .  jj./l  .176655  )*SQRT(  6.  37369+ (.  03*Y("i))**2)  *SQRT(1.  +,  0000027 


1 8 5T8~4~7Tl *  y  ( I )  * *2"j 

WRlTEC3,112)L»X/8,sSGM,sET 

GOT  0  3 0  0 

130  X00=> (I )/i .129333 

Sigma  =  c"i.  oT/i«  129333  >  *5qrt  ( 1 . 5  9  028 + ( . q  3  * y  <  m  *  *  2  >  *  SQfH  (1.  +  . 0000030' 

124294283*_MI)**2) 

w R"rf"r("37iT2  Ti»  x~o it,  sign  a  ,  set 

GOTO300 

150  X86=Y(I)/1. 182663 

SIGH  A  =  ( 1  .  (j  1  /  1  . 1  "2  663  )  •  SQR  T  ( 1  ,  9693  3+  (  .  03  +  Y  ( I J  )  *  *2  )  +  SQR  1  ( 1  .  +  .  Q  0  0  0  0  2  7_ 
157694084*r(I)**2) 

W  R  I T  E  (  3  ,  1 .1  2  )  L  »  X  8  6 ,  S I G  f 1  A  ,  S  E  T 

GOT 03 00 
36  0  X47sY(I)/l, 1433/0 

SIGMA=(l.Ul/l.l4337  0)*SORTU.00434  +  (,OL<  +  Y(I))*  +  2)*5QRl(l.  +  .OOon029 

150488846* r (I)**2) 
WRITE  (3 , 11 2 >  L .  X„47# SIGM Aj  SET 
'SW    CONTINUE 

4  q  n  c  0  M  T I  N  u  E  

G  O  T  01 0  0  0 

END 
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6.        Data  Interpretation   III 

The  meaning  of  the  terms   in  this  program  are  given  by  the   equation 
in  program  2. 

)I*ENSIQN    Y(lu) 
l"jj8    FORMAT  <2IJ) 

HO    FQRMAT  (  3Fin  .  6) 

112    FQCMAT(  /»  IPX*  1WX#  13, 3^    =     ,  -  10  ,6,5x,  7H+    OR    -     , F1U . 6, 5X , 9HPATA    SF- I     » 
113) 


11/    FORMAT<5Fl0.6,Fx5,9,FlO,6) 
C  INTERPRETATION    OF    DaTa    FROM    2    R    DOSIMETERS 

10r;n    R£AD<1»108)N0EXP 

IF<N0£XP.E0. 0)CALLEXIT 

1)0400  j  =  l,  N0EXP 

REAO(lVi0a)<M"PTS 

READ(X,108)L,SET 

REAQ(l,ll'))(Y(i;,[  =  l,\PTS) 

READ (1, 114 )X AVE #YINT» SLOPE *TOVB, BUNK »SIGMAD#SIGMA2 

DOSffOI'si.NPTS 

XgXAV£+(  Y(D-  flNT)/SLOPE 

SlGMA-;OVb*SL:HTlSIGl"iA2+(,a3*Y(I))**2)*SORT(l.+BU\1K+(YiI)-YifviT)**2/ 


1SIGMAD) 
WRI  rE(3»112)L»X, SIGMA, SET 
300    CONTINUE 

4Q0    C 0 N T IN o e  

GOTOioob 

END 


7.        Data  Interpretation  IV 


The  terms  here  have  the  same  meaning  as  elsewhere  except   for  SIGMAD, 

which  in  this  program  is  equal  to  SIGMA1  in  program  3  and  is  equal  to 

1 
SIGMAD 


in  program  2. 


DIMENSION    Y(iu) 
1  p  A    F  0  R  M  A  i  (  ^  [  ,5 ) 
110    FORMAT(5F10,6) 
112    FORMAT  </,i0X*iHX7;t3"#3H    =     »  F  i  r'  .  6  ,  5X  ,  /h  +    OR    -     ,  Fl  0  .  6  ,  !>X  .  QHPAt  A    ^K  I     » 

113) 
114    FoRhA'T'Tii  i  ,  r  1 1  .  b  ,  I-  i  1.  .  h. ,  F  J  8  .  15  ,  F 1  0  .  6  ) 
C  [NTERPRETATION    UF    P-ATa    FRO^    2    *    DOSIMETERS    i^OHlTO    Ti-»u    ORIGIN) 

100  0    Rkaq( 1» 108)NO£XP 


ioU 


IF  (NQEXP.EQ",  0  >CALLEXIT 
DQ4Q0j=l,NOEXP 
RE'AD"(lVl~n8)Nl»TS 
READ<l»108)LtS£T 
REA'D(iViiG)  ( Y(I)  ,  i=T»NPtS) 

re ad( i *ii*) slope* sigma2,sigm ad »tovb 

D030  0Isl,NPTS 

X  =  Y(  D/SL03E 

SIGMA  =  ]O^B*bQHTCSinr-iA^+(.03*Y(I))**2)+SURT(l.+SlGMAD+(Y(I)*  +  2)> 

WRITE(3*112)LiXj^SI6MA,SEJ 
30  0    CONTINUE 
4oO    CONTINUE  

Gcncnooo 

END 


8.   Data  Normalization 


NOEXP  The  total  number  of  sets  of  data  (3U) 

M  The  number  of  the  set  of  data 

NE  The  number  of  detector  locations 

NPTS  The  number  of  repetitions  of  the  experiment 

CURIES  Source  strength  in  curies 

SIGMAC  Error  term  associated  with  the  source  strength  ($>) 

D  Dose  rate  (mr/hr) 

S  Error  associated  with  the  dose  rate  (mr/hr) 

SIGMAD  Error  associated  with  the  dose  rate  ($) 

F  Normalization  factor 

DOS  Normalized  dose  rate  (mr/hr  per  curie/ft  ) 

SIG  Error  associated  with  the  normalized  dose  rate  (mr/hr  per  curie/ft  ) 

AVEDOS  Normalized,  averaged  dose  rate  for  a  particular  detector  position 

SIGMA1  The  normalized,  averaged  standard  deviation  associated  with  AVEDOS 
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DIMENSION    D<5)iSIQMAD(5)7t"lME"(5)#F(5)#D0S(5)»8IG(5y#AREA(5")»S(5") 
08    F0RMAT/(25X,8hAVED0S  =    »Fii.6) 
99    F0RMAT(/*25X/7HSIGMA=    » F 1 0  .  6  #  //"j 

100  f 0  r m a  r ( 1 j  ) 

1 0 1  F  0  R  M  A  T  ( 5  F  i  0  .  6 ") 

1  n  T5    FORMAT  (i0X#6Fi2.5#/) 

1  r-  4    F  0  R  M  A  T  ( 5  F 1 2  ,  5  j 

105    F0RMAT(5X»9HDATA    SEl     ,16) 

DAlA    NORMALIZATION     TO    R/HR/TcWlC/F  Y2) 

READ<1*100)NOEXH 

ifTnoexpteq.  orcAii  EXIT 

DO  230  j*l,NOEXP 


RLADd.  100  >M 
READ<  j.  '  10¥)NE 

read<  1  #  1 0 1> )  ,m  p  rs 

RE  AD  ( 1  *  1 0  1 )  OTMTC I )  *  1  =  1  f  NPTS ) 

R  E  A  0  (  1  ,  1  H  )  (  A  R  F.  A  (  1  )  .  1  =  1  ,  \  p  T  5  ) 

REaD  (  3  ,  10 X  JCu'WItsVslGMA'C 

PTSsNPTS 

Oq  208  M=3 #nE 

READ ( 1*101) <D(I)#I=l*NPTS> 

READ(  1,101)  (S(i)#i  =  i#NPfS"j 

00  220  K=i,\P,s 

F <K Y = A R'.E A '  <  K ) 7 T C U R IE S i * T I M E  (<)/ 60" ,") 


22'J 


1^0 


DOS(K) =F<K)*D(K) 
SIGMAD(K)=S(K)/D(K) 

SlG(k<)=DOS(K)«SuRT(SToMAn(<)*4.2 

I F ( D ( K ) . NE . 0 . )  b  0  TO  2  2 0 

S I R  (  K  )  =  0  . 

DOS(K)=0 , 

CONTINUE 

SUMDOSsQ . 

SIGMA 2= 0 . 

DO  180  L= 1  ,NP  1  S 

SIGMA2=SIGMA2  +  SIG(L)**2 

SuMDOS-SO<vriuS  +  DuS(l    ) 


♦     S  1 0  M  A  C  ♦  *  ?  ) 


2S0 


AVFDOS-SU-  nuS/P  I  S 

S I G  M  A 1  s  (  S  Gl  R  T  (  S I  li  M  A  2  )  V  /  P  t  S 

WRT  rt  (3,105)M 

MRlT£(4fli)3)(F(K)»K«l«NPTS) 

WRITE (3*1  3><DOS<K)#K=l»NPTS> 

WRITE(3,1u3;(SIIj(K)#Kb1#NPTS) 

WRIT*  (  i#98)AVLD0S 

WRITE  (3  ,;vv  )SIGMAJ 

WRITe(2#1U4)(d05(K)#Ks1#NI 

COM! I  JOE 

CONTINUE 
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ABSTRACT 

An  experimental  study  has  been  made  to  determine  the  radiation  intensity 
from  surrounding  fallout  at  various  positions  in  a  1:12  and  a  l:k   scale, 
steel  model  of  a  concrete  block  house  with  a  basement.   The  plane  fallout 
field  was  simulated  by  pumping  an  80-curie  cobalt-60  source  at  a  constant 
speed  through  plastic  tubing  laid  so  as  to  maintain  a  constant  tubing  length 
per  unit  surface  area.   Total  dose  measurements  were  made  using  ion  chambers. 
A  unique  feature  of  the  two  models  was  that  each  could  be  adapted  to  simulate 
a  block  house  with  a  portion  of  its  basement  walls  exposed. 

The  inverse  of  the  factor  by  which  the  radiation  intensity  is  reduced 
at  a  point  is  called  the  protection  factor  at  that  point.  With  the  aid  of 
a  method  developed  by  Kaplan,  et.  al.  the  experimental  data  were  analyzed 
and  protection  factors  were  calculated  for  locations  at  the  sides,  corners, 
and  center  of  each  model.  The  results  for  the  two  models  agreed  within  their 
standard  deviations  in  every  location  except  the  corners  of  the  basement.   A 
comparison  was  also  made  among  the  model  protection  factors,  those  calculated 
using  standard  computational  methods,  and  the  results  obtained  for  the  full- 
scale  building. 


